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This  dissertation  considers  issues  regarding  stochastic  systems  with  limited  re- 
sources. We  study  stochastic  multi-item  inventory  models  with  limited  warehouse 
capacity  for  cases  of  equal  and  unequal  replenishment  intervals.  We  also  discuss 
stochastic  production  systems  in  which  a single  machine  is  used  for  a series  of  cus- 
tomer orders. 

We  first  characterize  a stochastic  multi-item  inventory  model  for  the  case  of 
equal  replenishment  intervals  with  limited  warehouse  capacity  by  using  dynamic 
programming.  We  use  induction  to  derive  the  convexity  of  objective  cost  function 
and  show  that  a myopic  replenishment  policy  is  optimal  for  both  finite  and  infinite 
planning  horizons. 
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For  the  unequal  replenishment  interval  case,  we  completely  characterized  the 
optimal  replenishment  policy  for  the  two-item,  two-period  case.  For  more  general 
cases,  because  of  the  inability  to  show  convexity  of  the  objective  function,  we  con- 
struct three  heuristics  to  determine  replenishment  quantities  for  the  case  of  unequal 
replenishment  intervals.  Numerical  testing  of  these  heuristics  suggests  that  they 
yield  near  optimal  solutions  for  a small  set  of  problems.  The  performance  of  these 
heuristics  was  evaluated  for  a number  of  problems  with  various  cost  and  demand 
parameters. 

For  stochastic  production  systems  with  limited  resources,  we  have  modeled  the 
problem  of  setting  optimal  planned  order  leadtime  in  a single  stage,  make-to-order 
production  system.  As  a solution  approach,  we  derived  the  distribution  of  actual  pro- 
cess completion  time  of  individual  orders  and  then  determined  the  optimal  planned 
leadtime. 
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CHAPTER  1 
INTRODUCTION 


Supply  chain  management  has  become  an  important  principle  as  industries  strive 
for  operational  excellence  to  optimize  resources  and  to  reduce  costs.  By  running  an 
effective  supply  chain,  companies  can  meet  their  customer  demand  satisfactorily  and 
improve  their  competitive  position  in  markets.  One  of  the  tools  that  lends  itself  to 
improving  SCM  is  outsourcing.  That  is,  you  must  determine  when  to  have  some- 
one else  do  work  for  you  such  as  warehousing,  product  shipment,  or  manufacturing 
components,  so  that  you  can  concentrate  on  your  core  competencies  -those  things 
that  you  do  best.  For  example,  Carrier’s  McMinville  manufacturing  plant  (Swenson 
1998)  developed  a partnership  with  a supplier  of  non-production  service  parts.  This 
supplier  staffs  McMinville’s  service  parts  warehouse  and  is  fully  responsible  for  pur- 
chasing stocking,  and  scheduling  service  parts.  The  supplier  can  source  service  parts 
and  actually  run  the  warehouse  more  cost-  effectively  than  McMinville  can  with  its 
own  employees.  Federal  Express’s  warehousing  and  distribution  service  for  several 
manufacturers  (such  as  Dell  Computer,  Hewlett-Packard  and  National  Semiconduc- 
tor) is  another  outstanding  example  of  outsourcing.  Where  once  the  transportation 
process  was  seen  as  a lost  cause,  it  is  now  viewed  as  an  important  part  of  supply 
chain  where  business  is  focusing  on  reducing  the  cost  of  goods.  When  coupled  with 
warehousing,  the  distribution  process  commands  respect.  Federal  Express  is  cur- 
rently responsible  for  managing  the  aforementioned  companies’  products  warehouse 


1 


2 


as  well  as  the  entire  distribution  process  from  those  manufacturers  to  their  respective 
customers.  Continuous  improvement  in  reducing  average  delivery  time  by  companies 
such  as  Federal  Express  and  UPS  becomes  a critical  piece  of  cost  reduction  to  those 
manufacturers. 

On  the  other  hand,  on  the  service-providers’  side,  for  example,  Federal  Ex- 
press runs  a giant  warehouse  to  accommodate  those  products  from  outsourcing- 
partnerships  with  manufacturers  as  well  as,  of  course,  its  primary  business  such  as 
mail  and  parcels.  Traditionally,  warehouses  function  mainly  to  store  products  sup- 
plied from  manufacturers  or  suppliers.  Manufacturers  decide  when,  for  how  long,  and 
how  many  products  are  stored.  Warehouses  themselves  have  no  control  over  supply, 
but  mainly  concentrate  on  warehouse  or  storage  management.  Today,  instead  of 
functioning  merely  as  service  providers  of  products  warehousing  and  distribution, 
these  companies  could  play  a central  role  in  a supply  chain  in  which  they  are  in- 
volved. In  other  words,  once  the  service  providers  develop  a strong  relationship  with 
their  partners,  they  might  expand  their  roles  not  only  to  play  as  current  service- 
providers  but  also  to  control,  at  least  partially,  manufacturers’  production  schedules. 
This  could  be  done  by  offering  counter-compensation  to  manufacturers  through  a 
negotiation  process,  although  an  interesting  issue  - how  to  negotiate  - is  not  covered 
in  this  study.  For  example,  Federal  Express  might  offer  a discounted  transporta- 
tion cost  to  Dell  Computer  as  compensation.  Once  these  service  providers  are  in 
the  driver’s  seat  in  their  supply  chain,  they  can  contact  their  partners’  customers 
directly  and  collect  information  about  demand.  Based  on  this  customer  information, 
service  providers  can  ask  their  partners  to  produce  a certain  amount  of  product  to 
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meet  future  demand. 

Once  in  a position  to  play  a central  role  in  their  total  supply  chain,  and  as  the 
number  of  partnerships  grow,  service  providers  confront  how  to  run  the  entire  system 
effectively  with  their  own  limited  resources,  such  as  warehouses  and  transportation 
resources.  For  example,  Federal  Express  might  need  a well-defined  strategy  to  accom- 
modate products  from  its  partners  effectively  and  to  meet  their  customer  demands, 
given  that  its  warehouse  has  a limited  capacity  (see  Figure  1.1).  In  other  words, 
Federal  Express  may  need  to  determine  optimal  replenishment  policies  for  minimiz- 
ing costs.  Meanwhile,  the  company  maintains  partnerships  with  several  different 
manufacturers  at  the  same  time  and  each  of  them  likely  has  a distinct  production 
schedule.  The  company’s  warehouse  thus  would  need  a policy  of  reviewing  its  inven- 
tory irregularly  or  non-periodically  instead  of  running  a conventional  periodic-review 
policy. 

When  individual  manufacturers  have  the  same  replenishment  interval  for  prod- 
ucts, and  thus  the  system  has  a periodic-review  inventory  policy,  deriving  an  optimal 
replenishment  policy  would  not  be  difficult  (Chapter  2).  However,  when  each  manu- 
facturer has  different  replenishment  schedules  or  has  unequal  replenishment  intervals, 
the  periodic-review  policy  is  not  valid  anymore.  It  would  be  difficult  to  determine  an 
optimal  replenishment  policy.  In  this  study,  we  focus  on  unequal  replenishment  in- 
tervals (either  different  replenishment  schedules  or  unequal  lengths  of  replenishment 
intervals).  We  try  to  determine  the  optimal  replenishment  policy  for  minimizing 
warehouse  costs. 

Stochastic  production  and  inventory  models  involving  deterministic,  capac- 
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Manufacturers  Warehouse  Customers 

Figure  1.1:  A multi-product  inventory  system. 

ity  (production)  constrained  periodic-review  policies  have  attracted  the  attention 
of  many  researchers.  Evans  (1967)  first  considers  this  issue,  in  which  he  models 
periodic-review  production  and  inventory  systems  with  multiple  products,  random 
demands  and  a finite  planning  horizon.  He  develops  the  form  of  the  optimal  policy 
for  the  multi-product  control  appropriate  for  the  system.  Since  then,  much  of  the 
literature  has  studied  periodic-review,  single-product  systems  with  production  ca- 
pacity constraint.  Florian  and  Klein  (1971)  and  De  Kok  et  al.  (1984)  characterize 
the  structure  of  the  optimal  solution  to  a multi-period,  single-item  production  model 
with  a capacity  constraint.  Federgruen  and  Zipkin  (1986a, b)  show  that  a modified 
base-stock  policy  is  optimal  under  the  discounted  cost-criterion  and  the  average  cost- 
criterion,  respectively,  both  for  finite  and  infinite  planning  horizons.  The  modified 
base-stock  policy  is  such  that  when  initial  stock  is  below  a certain  critical  number, 
produce  enough  to  bring  total  stock  up  to  that  number,  or  as  close  to  it  as  possible, 
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given  the  limited  capacity;  otherwise,  do  not  produce.  They  also  characterize  the 
optimal  base-stock  levels  in  several  ways.  Kapuscinski  and  Tayur  (1998)  provide  sim- 
pler proof  of  optimality  than  Federgruen  and  Zipkin  (1986a)  for  the  infinite-horizon 
discounted  cost  case,  based  on  results  from  Bertsekas  (19988).  Ciarallo  et  al.  (1994) 
and  Wang  and  Gerchak  (1996)  analyze  a production  model  with  variable  capacity 
in  a similar  environment  to  that  described  by  Federgruen  and  Zipkin  (1986a).  Wang 
and  Gerchak  (1996)  also  incorporate  variable  capacity  explicitly  into  continuous  re- 
view models. 

DeCroix  and  Arreola-Risa  (1998)  study  an  infinite-horizon  version  of  the  multi- 
product case.  They  establish  the  optimal  policy  for  the  case  of  homogeneous  prod- 
ucts, and  propose  a heuristic  policy  for  heterogeneous  products  by  generalizing  the 
optimal  policy  for  the  homogeneous  product  case.  Glasserman  (1996)  addresses  a 
similar  problem  to  DeCroix  and  Arreola-Risa  (1998)  in  a continuous-review  system. 
He  presents  a procedure  for  choosing  base-stock  levels  and  capacity  allocations  that 
are  asymptotically  optimal,  but  assumes  that  a fixed  proportion  of  total  capacity  is 
dedicated  exclusively  to  each  product.  The  use  of  asymptotic  analysis  is  similar  in 
the  spirit  to  Anantharam  (1989).  His  static  allocation  problem  contrasts  with  the 
dynamic  scheduling  problem  addressed  in  Wein  (1989)  and  Zheng  and  Zipkin  (1990) 
and  the  priority  scheme  in  Carr  et  al.  (1996).  Lau  and  Lau  (1988,  1996)  present 
formulations  and  solution  procedures  for  handling  a multi-product  newsboy  problem 
under  one  resource  and  multiple  resource  constraint.  Nahmias  and  Schmidt  (1984) 
also  investigate  several  heuristics  for  a single-period,  multi-item  inventory  problem 


with  a resource  constraint. 
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A multi-product  dynamic  nonstationary  inventory  problem  with  limited  ware- 
house capacity  is  first  considered  in  Veinott  (1965).  He  provides  conditions  that 
ensure  that  the  base  stock  ordering  policy  is  optimal  in  a periodic-review  inventory 
system  with  a finite  horizon.  Ignall  and  Veinott  (1969)  show  in  the  stationary  case 
that  a myopic  ordering  policy  is  optimal  for  a sequence  of  periods  for  all  initial  in- 
ventory levels.  Beyer  et  al.  (2000a, b),  use  the  dynamic  programming  approach  to 
derive  the  optimality  of  the  ordering  policy  for  the  average  cost  problem  and  to  prove 
the  convexity  of  the  objective  cost  function  as  we  did  in  this  study  coincidentally. 
Beyer  et  al.  extensively  also  prove  the  optimality  of  the  modified  base-stock  policy 
in  the  discounted  cost  version  of  the  problem.  In  Federgruen  and  Zipkin  (1986a), 
they  consider  both  production  and  inventory  capacities  with  single-product  systems, 
and  derive  the  convexity  of  the  expected  cost  function. 

Anily  (1991)  and  Gallego  et  al.  (1996)  study  a multi-item  replenishment  problem 
with  deterministic  demand.  Anily  (1991)  investigates  the  worst-case  behavior  of  a 
heuristic  for  the  multi-item  replenishment  and  storage  problem  and  derives  a lower 
bound  on  the  optimal  average  cost  over  all  policies  that  follow  stationary  demand 
and  cost  parameters.  Gallego  et  al.  (1996)  consider  two  economic  order  quantity 
models  where  multiple  items  use  a common  resource:  the  tactical  and  strategic 
models.  They  derive  a lower  bound  on  the  peak  resource  usage  that  is  valid  for  any 
feasible  policy,  use  this  to  derive  lower  bounds  on  the  optimum  average  cost  for  both 
models,  and  show  that  simple  heuristics  for  either  model  have  bounded  worst-case 
performance  ratios.  Other  literature,  e.g.,  Rosenblatt  and  Rothblum  (1990),  Goyal 
(1978),  Hartley  and  Thomas  (1982),  Jones  and  Inman  (1989)  and  Dobson  (1987) 
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studied  deterministic  inventory  models  with  warehousing  constraint. 

Although  much  literature  is  devoted  to  multi-item,  periodic-review  systems  with 
a production  capacity  constraint,  little  has  been  done  for  stochastic  inventory  models 
with  a warehouse-capacity  constraint.  In  particular,  to  our  knowledge,  no  one  has 
considered  the  case  in  which  products  from  different  manufacturers  have  distinct 
replenishment  schedules  or  unequal  lengths  of  replenishment  intervals  (we  call  this 
the  case  of  unequal  replenishment  intervals).  Unfortunately,  except  for  very  small 
problems,  it  is  difficult  to  obtain  an  optimal  replenishment  policy  in  this  case  because 
we  are  unsure  whether  the  objective  cost  function  is  convex.  In  this  study,  we 
develop  three  heuristics  to  estimate  replenishment  quantities  for  the  case  of  unequal 
replenishment  intervals. 

We  first  characterize  a stochastic  multi-item  inventory  model  for  the  case  of 
equal  replenishment  intervals  with  limited  warehouse  capacity  by  using  dynamic 
programming.  We  use  induction  to  derive  the  convexity  of  objective  cost  function  and 
determine  the  optimal  replenishment  policy  for  both  a finite  and  an  infinite  planning 
horizon.  From  this  derivation,  we  realize  that  when  the  warehousing  constraint  is 
tight,  the  overall  warehouse  capacity  is  divided  into  several  terms  each  of  which  is 
related  to  products  from  an  individual  manufacturer,  respectively,  and  is  a constant. 
We  consider  each  term  as  the  capacity  of  the  corresponding  product.  This  result 
leads  us  to  propose  heuristics  to  determine  replenishment  policies  for  the  case  of 
unequal  replenishment  intervals. 

In  the  first  heuristic,  we  determine  the  capacity  of  each  item  separately  and 
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apply  this  over  the  entire  planning  horizon.  Given  the  separate  capacity,  this  multi- 
item, unequal  replenishment  interval  inventory  model  is  reduced  to  multiple  single 
item,  equal  replenishment  interval  problems  and  we  are  able  to  derive  the  optimal 
replenishment  quantity  of  each  item.  The  second  heuristic  is  similar  to  the  first 
one  in  the  sense  that  we  use  a separate  capacity  for  each  item.  Instead  of  using 
a fixed  separate  capacity  for  each  item  over  the  planning  horizon,  this  heuristic 
adjusts  individual  capacities  in  each  period  according  to  the  availability  of  any  unused 
capacities  of  other  items.  Thus,  compared  to  the  first  heuristic,  it  provides  more 
room  to  those  items  that  are  to  be  replenished  and  we  expect  that  this  heuristic 
has  better  performance  than  the  first  one  since  it  has  more  opportunity  to  replenish 
products.  In  the  last  heuristic,  we  use  the  total  available  capacity  of  the  system  at 
the  moment  to  replenish  products.  Thus,  this  heuristic  always  maintains  the  greatest 
capacity  for  replenishment  among  the  three  proposed  heuristics.  We  then  compare 
the  performance  of  each  heuristic  using  numerical  examples. 

In  supply  chain  management,  limited  resources  in  individual  supply  chain  mem- 
bers critically  affect  modeling  of  the  supply  chain.  So  far,  we  studied  a stochastic 
inventory  model  with  a limited  resource  of  warehousing  capacity.  We  next  discuss 
another  stochastic  system  with  limited  resources.  Consider  a stochastic  production 
system  in  which  a single  machine  serves  a series  of  customer  orders  to  produce  a single 
item.  Because  of  a single  production  line  coupled  with  a single  machine,  the  amount 
of  production  in  a given  time  is  limited.  Customer  orders  arrive  in  this  single-stage, 
make-to-order  production  system,  through  which  individual  orders  need  to  be 
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processed  with  their  own  due-date.  Most  of  the  due-date-setting  models  assume  that 
the  due  dates  for  individual  orders  are  set  entirely  exogenously.  However,  in  certain 
practical  contexts,  each  order  arrives  with  a due  date,  indicating  some  future  time 
when  the  customer  wishes  to  receive  the  goods  ordered  and,  in  most  cases,  due  date 
setting  is  negotiable  and  is  the  responsibility  of  the  marketing  personnel  of  the  firm. 
For  the  marketing  group,  who  have  knowledge  of  customer’s  wishes,  it  is  important 
to  know  their  manufacturing  personnel’s  perspective  on  the  order  before  negotiating 
the  due  date  with  their  customer.  That  is,  they  would  like  to  know  the  optimal 
customer  order  leadtime  to  quote  for  the  order.  By  customer  order  leadtimes,  we 
mean  the  time  from  a customer’s  order  until  the  due  date. 

Uncertainties  in  manufacturing  environments  play  important  roles  in  determin- 
ing planned  customer  order  leadtimes.  These  uncertainties  are  mostly  due  to  com- 
plicated production  processes,  random  yields,  and  high  quality  requirements.  In- 
corporation of  these  factors  into  the  selection  of  planned  leadtimes  represents  an 
important  step  toward  increasing  the  robustness  of  manufacturing  planning.  The 
leadtime  uncertainty  itself  affects  many  aspects  of  costs  and  control.  It  is  particu- 
larly problematic  in  queueing  systems  because  tardiness  in  earlier  orders  may  delay 
subsequent  orders. 

The  flood  of  literature  in  this  area  focuses  on  leadtime  uncertainties.  Whybark 
and  William  (1976)  compared  safety  stock  and  safety  time  as  alternatives  for  buffer- 
ing and  found  that  safety  time  coped  better  with  leadtime  uncertainty.  Melnyk  and 
Piper  (1985)  examined  the  impact  of  PLT  magnitude  on  manufacturing  performance 
and  found  that  larger  PLTs  improved  on-time  deliveries.  St.  John  (1985)  found  that 


10 


larger  PLTs  increased  total  costs  and  supported  the  view  that  PLTs  must  be  as 
small  as  possible.  In  particular,  St.  John  showed  that  the  carrying  costapplied  to 
raw  materials,  work  in  process,  and  finished  goods  inventorycontributed  the  greatest 
amount  to  the  increase  of  the  total  costs.  In  contrast,  Marlin  (1986)  suggested  that 
inventory  levels  increased  with  the  leadtime  error  magnitude  rather  than  with  the 
PLT  magnitude.  Marlin  explained  that  this  contradiction  is  due  to  shortcomings 
in  the  simulator  used  in  St.  John’s  study.  Other  literature  such  as  Weeks  (1981), 
Kanet  (1986),  Yano  (1987a,  b,  c)  and  Moran  (1998)  studied  empirical  investigations 
of  PLTs.  Weeks  proposed  a one-period  inventory  model  to  set  the  optimal  PLT  and 
suggested  that  the  PLT  is  a function  of  the  cost  of  tardiness,  earliness,  due  date 
length,  and  the  manufacturing  lead-time  distribution.  Weeks  also  showed  that  the 
single  stage  model  with  tardiness  costs  is  equivalent  to  a simple  newsboy  problem. 

Some  efforts  have  focused  on  the  impact  of  the  leadtimes  and  their  variability 
on  optimal  inventory  decisions  and  system  performance.  The  reduction  of  lead- 
times  and  their  variability  is  a key  element  of  process  improvement,  which  in  turn 
is  at  the  core  of  the  Just-In-Time  approach  (Schonberger  1982,  1986  and  Shingo 
1989).  Porteus  (1985,  1986)  and  Zangwill  (1987)  examined  the  effects  of  setup-cost 
reductions  in  deterministic,  EOQ  type  production  systems.  Song  and  Zipkin  (1996) 
investigated  the  joint  effect  of  leadtime  variance  and  lot  sizes  in  a parallel  process- 
ing environment.  Hariharan  and  Zipkin  (1995)  extended  basic  inventory  models  to 
allow  random  demand  leadtimes  as  well  as  random  supply  leadtimes.  On  the  other 
hand,  Duenyas  (1995)  formulated  the  problem  of  quoting  optimal  leadtimes  as  a 
semi-Markov  decision  process.  He  considered  decision-making  on  leadtimes  from  the 
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customers’  perspective  instead  of  the  common  approach  of  considering  the  manufac- 
turers’ preference  (Duenyas  and  Hopp  1995). 

This  study  is  stimulated  by  Yano’s  work  on  PLTs  in  serial  production  systems 
(1987a).  Yano  developed  a solution  approach  to  solve  the  two-stage  problem  and 
applied  it  to  iV-stage  serial  systems.  The  limitation  of  their  approach  was  that  they 
did  not  address  the  effects  of  queueing  or  capacity  explicitly  in  the  model.  For  in- 
stance, Yano  assumes  that  an  order  that  has  completed  processing  in  one  stage  may 
be  held  until  its  due  date  (i.e.  no  early  delivery  to  customers  is  allowed).  However, 
it  is  also  assumed  that  once  that  order  arrives  at  the  next  stage,  it  is  processed 
immediately  without  being  delayed  by  preceded  orders.  This  assumption  greatly 
simplifies  their  solution  approach  for  determining  PLTs  in  serial  production  systems. 
The  assumption  used  by  Yano  is  not  pragmatic  in  many  situations.  Literature  such 
as  Duenyas  (1995)  and  Duenyas  and  Hopp  (1995)  considered  queueing  effects  in  de- 
termining PLTs.  In  this  study,  we  consider  the  aforementioned  delay  as  well  as  the 
possible  holding  cost  incurred  due  to  the  prohibition  of  early  delivery  to  customers 
in  single-stage  systems. 

This  study  addresses  single-stage  make-to-order  (MTO)  production  systems.  Our 
main  objective  is  to  determine  the  optimal  planned  customer  leadtimes  while  mini- 
mizing the  sum  of  expected  inventory-holding  costs  and  expected  penalty  costs  for 
exceeding  the  quoted  due  date,  so  as  to  quote  planned  leadtimes  to  customers  at 
the  time  of  order  arrival.  Since  the  area  of  modeling  leadtime  setting  where  demand 
is  sensitive  to  quoted  leadtimes  is  relatively  unexposed,  it  would  not  be  a critical 
shortcoming  to  assume  that  the  leadtimes  quoted  to  customers  are  independent  of 
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each  other.  Within  this  scope,  we  focus  on  modeling  the  appropriate  expected  costs 
under  a variety  modeling  assumptions  and  characterizing  the  optimal  policies. 

Our  approach  to  solving  the  problem  is  to  derive  the  distribution  of  actual  com- 
pletion time  of  the  process  for  an  individual  order  and  to  compare  this  to  the  cor- 
responding quoted  due  dates  to  obtain  the  expected  total  costs.  We  then  minimize 
costs  with  respect  to  the  decision  variable(s),  which,  in  this  case,  are  PLT(s).  We 
show  that,  in  special  cases,  the  single  stage  model  is  equivalent  to  a simple  newsboy 
problem  as  Weeks  (1981)  pointed  out. 

Next  is  an  outline  of  this  dissertation.  In  Chapter  2,  we  introduce  a stochastic 
multi-item  inventory  model  for  equal  replenishment  intervals  with  a limited  ware- 
house capacity.  We  use  induction  to  show  the  convexity  of  objective  cost  function 
and  then  determine  the  optimal  replenishment  policy  in  both  a finite  and  an  infinite 
planning  horizon.  In  Chapter  3,  we  relax  the  assumption  of  identical  replenishment 
schedules  for  items  and  characterize  the  corresponding  inventory  model.  We  derive 
the  convexity  of  objective  cost  function  and  determine  the  optimal  replenishment  pol- 
icy for  the  two-item,  two-period  problem  as  a special  case.  For  large  sized  problems, 
we  propose  three  heuristics  to  approximate  the  optimal  replenishment  quantities 
and  compare  the  performance  of  each  heuristic.  Chapter  4 addresses  the  problem 
determining  optimal  planned  leadtimes  in  single  stage  make-to-order  production  sys- 
tems. We  introduce  an  unique  solution  approach  by  deriving  the  distribution  of  the 
actual  process  completion  time  of  newly  arrived  orders.  In  Chapter  5,  we  discuss 
conclusions  from  our  study  and  briefly  summarize  future  research  directions. 


CHAPTER  2 

EQUAL  REPLENISHMENT  INTERVALS 
2.1  Introduction 

In  this  chapter,  we  study  a stochastic  multi-item  inventory  model  with  a warehouse- 
capacity  constraint.  The  demand  in  each  period  of  equal  length  is  a random  vari- 
able and  inventories  are  reviewed  periodically.  Veinott  (1965)  first  considered  this 
dynamic  nonstationary  multi-item  inventory  model  with  discounted  future  costs. 
Veinott  provided  conditions  that  ensure  that  the  base  stock  ordering  policy  is  opti- 
mal and  that  the  base  stock  levels  in  each  period  are  easy  to  calculate.  Ignall  and 
Veinott  (1969)  show  in  stationary  case  that  a myopic  ordering  policy  is  optimal  for 
a sequence  of  periods  for  all  initial  inventory  levels.  Beyer  et  al.  (2000a, b)  take  up  a 
concrete  version  of  the  abstract  model  considered  by  Veinott  (1965)  and  Ignall  and 
Veinott  (1969).  Beyer  et  al.  use  the  dynamic  programming  approach  to  derive  the 
optimality  of  the  ordering  policy  for  the  average  cost  problem  and  prove  the  convex- 
ity of  the  objective  cost  function.  They  present  a detailed  study  of  finite  and  infinite 
horizon  multi-product,  discounted-  and  average-cost  problems.  A part  of  their  work 
is  coincidentally  indistinguishable  from  what  we  study  in  this  chapter. 

We  first  characterize  a stochastic  multi-item  inventory  model  for  the  case  of  equal 
replenishment  intervals  with  limited  warehouse  capacity  by  using  dynamic  program- 
ming. We  then  use  induction  to  derive  the  convexity  of  objective  cost  function  and 
determine  the  optimal  replenishment  policy  for  both  a finite  and  an  infinite  planning 


13 


14 


horizon. 


2.2  Model  Assumptions  and  Formulation 
We  first  consider  single-item  inventory  systems,  facing  limited  warehouse  capac- 
ity. A warehouse  supplies  single  goods  from  a manufacturer  to  a population  of 
customers.  Customer  demand  in  each  period  is  an  independently  distributed  ran- 
dom variable.  When  the  demand  cannot  be  met  by  inventory  on  hand,  the  order  is 
backlogged  until  sufficient  inventory  is  available.  Distribution  of  the  demand  is  as- 
sumed to  be  stationary.  We  assume  that  the  leadtime  between  the  manufacturer  and 
the  warehouse  is  deterministic  and  is  set  here  as  zero  so  that  replenishment  orders 
from  the  warehouse  and  corresponding  shipments  occur  simultaneously.  Further,  a 
fixed  schedule  of  shipments  with  equal  intervals  is  preset  at  the  beginning  of  each  pe- 
riod. Therefore,  the  inventory  system  considered  here  is  identical  to  periodic-review 
systems.  As  mentioned  before,  we  do  not  assume  any  specific  inventory  policy  such 
as  the  base  stock  policy.  We  use  the  following  notation: 

Parameters 

m : number  of  replenishment  periods, 
j : index  of  replenishment  periods,  j = 1, 

V0  : fixed  storage  capacity  of  the  warehouse, 

T : length  of  the  time  horizon, 

h : unit  cost  of  holding  inventory  at  the  warehouse, 

p : unit  penalty  cost  for  demand  backordered  at  the  warehouse. 

Demands 

Dj  : demand  in  period  j , a random  variable, 
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State  and  decision  variables 

Ij  : inventory  level  at  the  beginning  of  period  j before  replenishment,  where 
1 1 is  given, 

Qj  : size  of  replenishment  in  the  jth  period, 


We  define  gj(Ij,Qj)  as  the  one-period  expected  holding  and  penalty  costs  in 
period  j,  given  the  initial  inventory  of  Ij  and  the  size  of  replenishment  of  Qj.  Then 
we  have 


9j(Iji  Qj) 


E[h(Ij  + Qi  - DjV  + p(Dj  - Ij  - Q,)+] 

roo 

/ [h(Ij  + Qj  - x)+  + p(x  - Ij  - Qj)+]  fj(x)dx 
J i=0 

rh+Qj  r°° 

h / (Ij  + Qj-  x)fj(x)dx  +p  (x  Ij  Qj)fj(x)d: 

Jx= 0 Jx=Ij+Qj 


(2.1) 


where  fj  is  the  probability  function  of  Dj. 

To  formulate  this  multi-period  inventory  model  by  using  stochastic  dynamic  pro- 
gramming (SDP),  we  define  Ij  as  the  state  variable  of  the  system  in  period  j.  The 
decision  in  each  period  is  to  choose  the  size  of  replenishment,  Qj,  for  periods  j 
through  m,  given  a starting  inventory  of  Ij.  Given  the  initial  state,  the  actions  and 
the  demands  in  period  j,  the  state  in  the  next  period  is  determined  recursively  as 
follows: 

Ij+i  = + Qj  ~ Djt  j — 1)  m (2-2) 

With  the  definitions  above  we  can  state  the  recursive  equations  whose  solution 
provides  the  optimal  policy  for  the  problem.  Let  Gj(Ij,  Qj)  be  the  expected  holding 
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and  penalty  cost  for  periods  j through  m along  with  initial  inventory  of  Ij.  Also, 

define  Gj(Ij)  = min  Gj(Ij,Qj ) as  the  minimum  expected  costs  in  period  j 
Qj  -Ij +Qj  < Vo 

through  m,  given  that  the  system  starts  period  j in  state  Ij  and  is  subject  to  the 
capacity  constraint  Ij  + Qj  < V0.  Thus,  the  following  expression  gives  the  expected 
cost,  Gj(Ij,Qj),  in  period  j through  m: 

Gj(Ij,  Qj)  = 9j{Ij,  Qj)  + E [G*+1(ij+i)] 

poo 

= Qj)  + / [g;+1(/,-  + Qj-  £>,)]  HDjdDj 

Jo 

poo 

= 9j(Ij,Qj)+  [Gj+l{Ij  + q - D,)]  f{Dt)dD,  (2.3) 
J 0 

where  Qj  is  the  optimal  solution  in  period  j.  The  optimal  cost  for  m periods  then 
becomes  G\(I i)  and  can  be  calculated  recursively  in  backward  by,  for  j = 2, 

poo 

G-(Ij)  = a ,imn  Qj)  + / [G’+ M + Qj-  £>,)]  f(Dj)dDt 

Note  that  the  capacity  constraint  is  a linear  function  of  Qj.  Thus,  it  is  enough 
to  show  convexity  of  Gj(Ij,Qj)  in  Qj  in  order  to  obtain  the  optimal  solution  of 
this  constrained  optimization  problem.  As  we  mentioned  previously,  we  first  show 
convexity  of  Gm(Im,Qm)  in  the  final  period,  determine  the  corresponding  optimal 
solution,  and  then  plug  it  in  the  expected  cost  function  in  the  previous  period  (m—  1). 
We  repeat  this  process  until  the  beginning  of  the  time  horizon.  Before  we  derive 
convexity  of  the  m-period  cost  function,  we  first  solve  a single-period  problem. 

For  simplicity,  we  drop  the  subscript  j in  a single-period  and  rewrite  the  problem 


as  follows: 
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min  G(I,Q)  = g(I,Q) 


subject  to  I + Q <Vq 


Let  L = g(I,  Q ) + /x(7  + Q -Vq)  where  fi  > 0.  Applying  the  well  known  Karush- 
Kuhn-Tucker  (KKT)  conditions,  we  obtain  the  following  system  of  nonlinear  alge- 
braic equations  and  inequalities: 


(b)  7 + Q < Vo, 

(c)  n{I  + Q-  Vq)  = 0, 

(d)  n > 0, 

Before  we  solve  these  equations,  we  introduce  the  rule  for  differentiating  integrals 
called  Leibnitz’s  Rule  (Casella  and  Berger  (1990)). 

LEIBNITZ’S  RULE  : If  f(x,6),a(9)  and  b(9)  are  differentiable  with  respect  to  9, 


then 


Notice  that  if  a(9)  and  b(9)  are  constants,  we  have  a special  case  of  Leibnitz’s  Rule : 


Any  interchangeability  of  integral  and  partial  derivative  in  this  study  follows  this 


rule.  From  (2.1)  and  (a),  we  have 
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8L 

dQ 


£"(i,Q)  + £w  + Q-v, »)] 

d_  r n+Q 
dQ 


rl+Q  f oo 

h (I + Q-  D)f{D)dD  + p ( D-I - Q)f(D)dD 

. Jo  Jl+Q 


+ QQ  [^(J  +Q  - VO)} 

rl+Q  r oo 

h / f(D)dD-p  / f(D)dD  + p 
Jo  Jl+Q 

hF(I  + Q)  - p(l  - F(I  + Q))+  p 


(h  + p)F(I  + Q ) — p + /i  — 0 


(2.4) 


where  F(-)  is  the  cumulative  distribution  function  of  D.  By  taking  the  second 
derivative  on  L,  we  get,  from  (2.4), 

ff  T fi2  r)2 

w = wgiI’Q)  + WWI  + Q-v°)] 

= (h+p)f(I  + Q)  (2.5) 


Since  the  right  hand  side  in  (2.5)  is  always  nonnegative,  the  expected  cost  for  a 
single-period  is  convex  in  the  replenishment  size,  Q.  The  following  two  propositions 
are  immediate  from  above. 


Proposition  2.1  The  expected  holding  and  penalty  cost  for  a single-period,  G(I,  Q), 
with  a limited  capacity  is  a convex  function  of  the  replenishment  size,  Q. 


Proposition  2.2  Q*  = min 


F 1 (fcfe) 


Vh 


— I is  the  optimal  replenishment  policy 


in  a single-period  problem 
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Proof : Solving  (2.4)  and  (b)-(d),  we  get  Q*  (See  Appendix  A for  details).  □ 
Now,  we  show  convexity  of  Gj(Ij,Qj ) in  Qj.  Assuming  Gj(Ij,Qj ) is  twice  differ- 
entiable, we  can  then  derive  the  following: 


Proposition  2.3  The  expected  holding  and  penalty  cost  for  periods  j through  m, 
Gj(Ij,  Qj),  is  convex  in  the  replenishment  size  in  period  j,  Qj. 

Proof  : First,  we  take  the  first  and  second  partial  derivatives  of  Gj(Ij,Qj ) with 
respect  to  Qj  from  (2.3): 


dGj(Ij,Qj) 

dQj 

dQj2 


d-^+rA[G^{ii+i)]Mx^  (2-6) 


The  first  term  in  the  right-hand  side  in  (2.7)  is  equivalent  to  the  second  derivative  of 
expected  cost  for  a single-period  problem,  and  thus  it  is  always  nonnegative,  referring 
to  Proposition  2.1. 

Now,  we  show,  using  induction,  that  the  second  term  in  (2.7)  is  nonnegative. 
First,  consider  period  m.  Note  that  this  SDP  problem  in  period  m is  identical  to  the 
single  period  problem  described  previously  because  G*m+l(Im+ 1)  = 0.  That  is,  the 
decision  in  the  final  period  m is  identical  to  that  for  a single  period  problem  since  the 
recursive  part  becomes  zero.  Since  the  first  and  second  derivatives  of  G^+1(/m+i) 
are  also  zero,  it  is  obvious,  from  Proposition  2.1,  that  Gm(/m,Qm)  is  convex  in  Qm. 
In  addition,  we  can  determine  the  optimal  solution  in  period  m,  based  on  the  result 
from  the  Proposition  2.2,  as  follows: 
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Qrn  = min 


7-1 


P 


h + p 


,V0 


~Ir, 


(2.8) 


Next,  we  assume  that,  for  period  (k  + 1),  Gk+i(h+i,  Qk+i)  is  convex  in  Qk+i  and 
the  optimal  solution  in  period  ( k + 1)  is 


Qk+i  = min 


p-i  ( JP_ 

k+1  ' h + p 


,V0 


-h 


k+l 


Substituting  Q*k+l  into  the  expected  cost  in  period  (k  + 1),  we  get 


jfc+i 

E [G*k+l(Ik+i)]  = E 
= E 


h (h+  i + Ql+i  ~ Ek+i)  + +P  (Dk+ 1 - h+ i - Qk+i) 


h min 


( t+z  i .Vo 


h + p J 

+ p ( Dk+ 1 - min 


— Dk+i 


,V0 


+ 


We  can  notice  that  the  optimal  expected  cost  in  period  (k  + 1)  is  not  a function  of 
decision  variables  any  more.  That  is,  E [G£+1(/fc+i)]  is  a constant  over  any  Qj,j  = 
1,  ■••>  m. 

We  now  show  convexity  of  Gk(h , Qk)  in  Qk-  Taking  the  first  and  second  partial 
derivatives  on  the  expected  cost  in  period  k, 


dGk(Ik,Qk) 

dQk 


[GUi(Ik+] 


d 


9k(Iki  Qk) 


dQk ' 

(h  + p)Fk(Ik  + Qk)  ~ V 


d 


e [c:+1(/fc+1)]  =0 


dQk 

(from  (2.1)  and  (2.4)) 


d2Gk(Ik,Qk) 
dQk 2 


= (h  + p)fk(h  + Qk) 
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and  the  right  hand  side  in  the  second  derivative  is  always  nonnegative.  Therefore, 
Gk(h,Qk ) is  a convex  function  of  Qk ■ Finally,  by  induction,  we  conclude  that,  for 
any  period  j,  Gj(Ij,  Qj)  is  convex  in  Qj.  □ 


The  result  in  Proposition  2.3  shows  that  the  solution  satisfying  the  first-order 
condition,  (2.6),  minimizes  the  expected  holding  and  penalty  costs  for  periods  j 
through  m , Gj(Ij,Qj),  and  we  can  easily  obtain  the  optimal  solution  in  period  j, 
analogous  to  (2.8): 


Proposition  2.4  Q*  = min  F-  1 , V0  - Ij  is  the  optimal  replenishment  pol- 

icy in  period  j,  where  Ij  is  the  initial  inventory  at  the  beginning  of  period  j , j = 


1,  m. 


Alternatively,  we  can  express  the  optimal  solution  above  as  follows:  Q*  = 

Ffl[(p  - p*)/{h  +p)]  - Ij,  where  p*  - arg  min  Ffl[(p-  p)/(h  + p)]  < V0.  Ac- 
cording to  Proposition  2.4,  the  optimal  replenishment  quantity  in  one  period  is  a 
function  of  its  initial  inventory  as  well  as  the  distribution  function  of  its  demand  and 
the  storage  capacity.  Furthermore,  we  can  see  that,  referring  to  (2.3)  and  Propo- 
sition 2.4,  the  recursive  term  in  the  expected  cost  is  a constant  over  the  current 
decision  variable,  Qj . That  is,  the  future  cost  incurred  in  that  period  is  not  affected 
by  the  replenishment  quantity  determined  at  the  beginning  of  the  period.  This  result 
gives  us  a good  intuition  when  we  extend  this  finite  horizon,  periodic-review  model 
to  the  infinite  horizon  case  because  deriving  convexity  of  the  expected  cost  function 
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in  any  period  is  not  affected  by  how  many  replenishment  periods  are  remained  in 
the  future.  We  will  discuss  this  issue  in  detail  later. 

2.3  Multi-Product  Inventory  Systems 

In  this  section,  we  extend  the  single-item  inventory  system  to  the  multi-item  case. 
Consider  a multi-item  inventory  system,  in  which  a warehouse  serves  manufacturers 
by  storing  their  products  and  delivering  them  to  their  respective  customers.  Each 
manufacturer  produces  a single  item  which  is  shipped  to  the  warehouse  according 
to  its  preset  schedule.  The  replenishment  leadtimes  between  the  warehouse  and 
manufacturers  are  assumed  to  be  zero  so  that  the  replenishment  orders  from  the 
warehouse  to  manufacturers  and  the  corresponding  shipments  occur  simultaneously, 
and  all  replenishment  are  received  at  the  beginning  of  the  following  replenishment 
period.  This  situation  might  occur  in  a JIT  environment  when  suppliers  are  located 
in  close  proximity  to  a warehouse. 

Again,  we  assume  that  items  are  replenished  periodically  by  manufacturers.  Cus- 
tomer demand  for  items  in  each  period  is  an  independently  distributed  random  vari- 
able, and  is  backlogged  when  it  cannot  be  met  by  inventory  on  hand  until  sufficient 
inventory  is  available.  The  distribution  of  the  demand  is  assumed  to  be  stationary. 
Most  of  the  notation  used  here  is  analogous  to  that  in  the  single  item  case  except 
for  the  extra  index  for  items  as  follows: 

Parameters 


l : number  of  product, 
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hi  : unit  cost  of  holding  inventory  for  item  i, 

Pi  : unit  penalty  cost  of  item  i for  demand  backordered  at  the  warehouse. 
Demands 

Ditj  : demand  of  item  i in  period  j,  a random  variable, 

State  and  decision  variables 

Iij  : inventory  level  of  item  i at  the  beginning  of  period  j before  replenishment, 
where  is  given, 

Qij  : size  of  the  replenishment  of  item  i in  the  jth  period,  i = 1, ...,/, 

We  first  define  Sj  = (Iij,  •••,  Iij)  as  a vector  of  state  variables  of  the  system  in 
period  j and  let  Qj  = {Qij,  --,Qi,j)  be  a vector  of  decision  variables.  For  a given 
state  Sj,  the  current  inventory  level  of  the  system  in  period  j is  Ij  = Ii,j>  which 
can  be  expressed  recursively  as  follows: 

i i 

Ij  = Ij- 1 + ^ ] Qi,j- 1 — 'y  ' Dij,  j = 2, ..,  TTl  (2-9) 

i=l  i— 1 

We  also  can  determine  the  current  inventory  levels  of  individual  items  recursively 
using  the  state  in  the  previous  period: 

Ii,j  = Iij-1  + Qi,j- 1 — Di jt  j — 2,  ...,  TTl  (2.10) 

We  define  gj(Sj,Qj ) as  the  expected  holding  and  penalty  costs  of  the  system  in 
period  j,  given  initial  inventories  of  Sj  = (hj, Iij)  and  replenishment  sizes  of 
Qj  — Then  we  have 

i 

Pj(Sj,  Qj)  — ^ E [hi(Iij  + Qij  - Ditj)+  + Pi(Ditj  — Iitj  — Qy)+]  (2.11) 

i=l 
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We  state  the  recursive  equations  whose  solution  provides  the  optimal  policy  for 
the  problem.  Denote  the  expected  holding  and  penalty  costs  for  periods  j through 
m by  Gj(Sj,Qj).  Then  Gj(Sj,  Qj)  can  be  expressed  as  follows: 

Gj(Sj,Qj)  = ®(Si,<M  + E[G’+1(Sj+ 1)] 

roc  roc 

= 9j{Sj,Qj)  + •••/  [G'j+l(Sj+l)]  flJ(xl)---fltj{xl)dxl---dxl 

Jo  Jo 

(2.12) 


where  /y  is  the  probability  density  function  of  a random  demand,  Ditj,  G*k(Sk)  = 
min  Gk{Sk,Qk ) is  the  minimum  expected  costs  in  period  k through 

m,  given  the  system  starts  period  k in  state  5*  = (h,k,  h,k)  and  is  subject 

to  the  capacity  constraint.  We  eventually  want  to  find  the  optimal  cost  for  a whole 
planning  horizon,  which  can  be  calculated  recursively  in  backward  by,  for 

j = 2,...,m, 


G’(St) 


min 

Qi,yZli=i(iij+Qij)<Vo 


[G*+i(Sj+i)]  /i,j(xi)  • • • fitj{xi)dx i • • • dxi 

As  in  the  single-item  case,  the  capacity  constraint  is  a linear  function  of  Qitj  and 
thus  it  is  enough  to  show  convexity  of  Gj(Sj , Qj)  to  get  the  optimal  solution  of  above 
SDP  problem.  Applying  the  same  induction  method  as  in  the  single-item  case,  we 
have  the  following  two  propositions. 
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Proposition  2.5  The  expected  holding  and  penalty  cost  for  periods  j through  m, 
Gj(Sj , Qj ),  is  convex  in  replenishment  sizes  in  period  j,  Qij  where  Qj  = 

(Q  ij>  Qij)- 

Proof  : Taking  the  first  and  second  partial  derivatives  of  (2.12)  yields  followings: 
For  any  i,  i'  = 1, n, 


dGj ( Sj , Qij, Qij) 


dQij 

Qij) 

dQij2 

> Qij) 

dQijdQvj 


d 


dQ 


~9j{Sji  Qij,  Qij) 


hj 

oo 


+ 


/ ’•••/ 


00  dG*+l 


dQk 


fij(xi)---fij(xi)dxi---dxi 


(2.13) 


d 2 


dQ 


>9j {Sj , Qij,  •••>  Qij) 


+ 


hj 

poo  poo 

J 0 «/0 


J + 1 


fij(xi)  ■■  -fij{xi)dxi  ■••dxi 


d 2 


dQijdQ 


9j{Sj,  Qij,  -,Qij) 


ij 


/»oo  /*oo 

lo  '"/» 


a2G 


j+i 


dQijdQfj 


~fij{xi)---fij(xi)dxl---dxl 


The  interchangeability  of  integral  and  partial  derivative  above  is  based  on  Leib- 
nitz’Rule.  We  use  Hessian  matrices  to  express  above  second  derivatives.  Define: 


d2Gj/dQitj2  •••  d2Gj/dQitjdQitj 


Hi  = 


d2Gj/dQiJdQi>j 


d2Gj/dQijdQ\j 


d2Gj/dQij2 
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Pgj/dQif 


&9j/dQi,jdQld 


d2  9 j!  dQi,jdQi'd 


d29j/dQi,jdQid 


d2gj/dQif 


d2Gj+i/dQid2 


d2G*j+l/dhd+ldQld 


&G)+JdQijdQVj 


d2G*j+l/dQlddQld 


d2G*j+l/dQld2 


Then  Hessian  matrices  defined  above  have  the  following  relationship,  based  on  (2.12): 


Note  that  Bj  is  equivalent  to  Hessian  of  a single-period  expected  cost,  and  is 
always  positive  semi-definite  in  any  period  j (see  Appendix  B for  details).  We  use 
induction  to  show  that  Hj  is  also  positive  semi-definite.  Consider  the  final  period  m. 
Assuming  G^+1(5m+i)  = 0,  it  is  obvious,  by  definition,  that  Em+ 1 is  a zero  matrix 
and  Hm  = Bm.  Therefore,  Hm  is  positive  semi-definite  because  Bm  is,  or  Gm(Sm,  Qm) 
is  a convex  function  in  Qi,m-  Furthermore,  we  can  determine  the  optimal  solutions 
in  period  m using  the  convexity  of  Gm  (see  Appendix  C for  details): 


where  /rm  is  a Lagrange  Multiplier,  FijTn  is  the  cumulative  distribution  function  of 
fi,  and 


Hj  — Bj  + 


•••/  [Ej+i\fi,j(xi)-“fi,j{xl)dxi---dxl  (2.14) 


0 Jo 
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We  assume  that  Hk+\  is  positive  semi-definite,  and  that  the  optimal  solution  in 
period  (k  + 1)  is 


/-)*  ri-l  ( P*  Pjs+l  i _ t 

1 — ri,k+ 1 l u i ->  / 


hi  + Pi 


where 


Pk+i  = ar9  min 

M*+i>0 


|JFi'l+iV  ) 

Then  the  optimal  expected  cost  in  period  ( k + 1)  can  be  obtained  by  substituting 

Qk+i  hy  Qk+ 1 m (■S'fc+i ) Qfc-t-i)- 

G^+1(Sfc+1)  = Gk+i(Sk+uQl+1) 

i 

= Y,E[h>(IiMl  + Q'iMi-D,Mi )+ 

i=l 

+Pi  {Difi+l  ~ Ii,k+l  ~ Qi,fc+i) 

= ' 


-i  A 


< Vo 


(2.15) 


»=i 


A‘  ( w)  - A-') 

+ Pi  ( Di,k+\  — 


-1  I Pi  - Pk+l 


hi  + pi 


(2.16) 


We  can  see  that  the  optimal  expected  cost,  G*k+l(Ik+ 1)  is  a constant  over  any 
decision  variables.  We  now  show  positive  semi-definiteness  of  Hk.  From  (2.12),  we 
get  the  following: 


d 


dQi,k 


Q 

Gk{Sk,  Ql,ki  ".,  Ql,k)  = gk(Sk,  Ql,k,  Ql,k) 


dQi„ 


i.k 


d 


dQitl 


i.k 


E [G*k+l(Sk+l)] 


d 


dQi„ 


i.k 


~9k{Sk,  Q\}ki  •••)  Ql,k ) 


(■■■ E [g;+1(s*+i)]  = o) 
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d2  d2 

~(*k(Sk,  Q\,k>  ■■■>  Ql,k)  — 7^  %9k(Fk,  Ql,kt  •■•)  Ql,k) 


dQlk 

a2  92 

Gk{SkiQ\,kT--,Ql,k)  — -KTk 9k{Sk,  Qitk,  ■■■,  Ql,k) 


dQi,kdQi',k 


9Qi,k9Qi',k‘ 


By  definition,  Hk  = Bk  That  is,  Hk  is  positive  semi-definiteness  because  Bk  is. 
Therefore,  we  can  conclude  that  Hj  is  convex  in  Qitj  for  any  period  j,  j = 1, 

□ 


The  result  in  Proposition  2.5  shows  that  the  solution  satisfying  the  first  order 
condition  is  a unique  optimal  solution  for  item  i in  period  j.  The  solution  approach 
in  here  is  analogous  to  that  in  the  single-item  case  since  the  expected  cost  is  a linear 
function  of  decision  variables  (see  (2.11)),  and  the  recursive  term  is  a constant  over 
those  variables.  The  following  proposition  states  the  optimal  replenishment  policy 
for  item  i in  period  j in  multi-item,  periodic-review  inventory  system: 


Proposition  2.6  The  optimal  replenishment  policy  for  item  i in  period  j is 


hi  + Pi  J 


-I, 


(2.17) 


where  is  a Lagrange  Multiplier,  Fid  is  the  cumulative  distribution  function  of  a 
random  demand  Ditj,  and 


arq  mm 
^■>o 


(2.18) 
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According  to  Proposition  2.6,  the  optimal  replenishment  quantity  of  an  individ- 
ual item  in  one  period  is  a function  of  its  initial  inventory  as  well  as  the  distribution 
function  of  its  demand  and  the  storage  capacity.  Furthermore,  Proposition  2.6  leads 
us  to  state  that,  in  any  period  j,  Q*j  + /y,  which  is  the  optimal  inventory  level  of 
item  i after  replenishment,  is  a constant  over  the  decision  variable  in  that  period, 
Qij.  Therefore,  as  we  have  seen  in  (2.15)  and  (2.16),  after  replacing  Q*j  + Ii}j  by 
FJ“1(-),  the  recursive  term  in  (2.12),  G*+l(Sj+i),  becomes  independent  of  Qij.  When 
we  derive  the  convexity  of  the  expected  cost  function  in  any  period  j,  Gj(Sj,Qj), 
we  then  only  need  to  show  the  convexity  of  the  single-period  expected  cost  function, 
gj(Sj,Qj ) because  the  recursive  term,  £’[G*+1(5j+1)],  is  a constant  over  Qy.  And 
deriving  the  convexity  of  gj(Sj,Qj ) is  usually  not  difficult,  according  to  Appendix 
C.  As  we  mentioned  in  the  single-item  case,  this  result  would  be  valuable  when  we 
extend  this  finite  horizon,  periodic-review  model  to  a infinite  horizon  problem.  Be- 
fore we  move  to  an  infinite  horizon  case,  we  state  the  following  theorem  based  on 
above  arguments. 

Theorem  2.1  The  optimal  replenishment  quantity  in  a single  period  is  independent 
of  the  future  expected  cost  incurred  in  that  period 

2.4  Numerical  Results 

In  this  section,  we  illustrate  numerical  solutions  of  the  optimal  replenishment 
policy  obtained  in  (2.17)  for  various  cost  parameters  and  capacities. 
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Exponentially  distributed  demand 

If  demand  of  item  i in  period  j is  exponentially  distributed  with  a mean  1/A ij, 
the  cumulative  distribution  function  of  the  demand  is  as  follows: 


= 1 — e 


and  the  inverse  of  Fitj(x ) is 


Fij(y)  = ~Xij  log{l-y) 


From  (2.14)  and  (2.15),  we  get 


Qij  — Ajj  log 

a*  = arq  min 
3 N>° 


(hi+  hj\  _ r 

\hi  + Pi  ) ' 


■ l 

XI  “A id  log 
_i= 1 


/ hi  Pj  \ 
\hi  + Pi  J 


Numerical  solutions  for  Qi,j  in  above  equations  can  be  obtained  from  the  following 
example. 


Example. 

Let  l = 2 and  j — 1.  Let  the  mean  demand  of  each  item  for  one  period  be 
A 1,1  = 100,  A2,i  = 120,  respectively,  and  the  capacity  limit  is  = 200.  Assume  that 
the  initial  inventory  levels  are  /i,i  = 10  and  /2ii  = 20,  and  the  holding  and  penalty 
costs  are  hi  = 5,  h2  = 5,pi  = 10,  and  p2  — 10.  Then  we  get  //•  = 1.04.  Thus,  the 
optimal  replenishment  quantity  of  each  item  is  <5i,i  = 81  and  Q2  i — 89. 
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1.  Effect  of  capacity  on  replenishment  quantity. 

Figures  2.1  and  2.2  represent  the  relationship  between  replenishment  quantities  of 
individual  items  and  warehouse  capacities  for  different  unit  costs.  Not  surprisingly, 
these  figures  indicate  that  the  replenishment  quantity  of  each  item  increases  up 
to  a certain  level  as  capacity  increases  and  then  becomes  flat  afterward  for  any 
combination  of  cost  parameters.  That  is,  at  a certain  level  of  capacity  the  capacity 
constraint  is  not  binding  anymore  to  the  problem. 

We  also  notice  that,  in  Figure  2.1,  item  1 has  more  replenishment  quantity  than 
item  2 when  the  capacity  decreases,  although  the  latter  has  a greater  mean  demand. 
This  is  because  when  the  size  of  capacity  is  small,  the  difference  in  the  optimal 
inventory  level  (after  replenishment)  of  each  item  would  not  exceed  the  difference  in 
the  initial  inventory  level  of  each  item,  so  that  item  2 is  replenished  less. 

2.  Lower  bound  of  capacity  in  each  period. 

The  replenishment  quantities  of  items  are  stabilized  at  a certain  level  of  capacity 
(250  in  Figure  2.1  and  210  in  Figures  2.2,  which  are  the  lower  bounds  of  the  capac- 
ity for  each  case).  By  specifying  input  parameters  and  running  a simulation,  the 
lower  bound  of  capacity  in  each  period  can  be  easily  determined.  It  is  advantageous 
to  warehouse  managements  who  want  to  utilize  their  storage  more  effectively  over 
periods. 


2.5  Infinite  Horizon  Case 

We  now  extend  to  a infinite  horizon  case.  While  we  label  periods  as  time  to  the 
end  of  the  planning  horizon  in  the  finite  horizon  problem,  it  is  typical  that  the  process 
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-Q(i,i) 

■Q(2,j) 


Figure  2.1:  Optimal  replenishment  quantities  versus  capacity  (hi  = 5,pi  = 10 ,h2  — 
5,p2  = 10) 
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-0(1  j) 
■Q(2,j) 


Figure  2.2:  Optimal  replenishment  quantities  versus  capacity  (/ii  = 5,pi  = 5 ,h2  = 
5,P2  = 10) 
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starts  at  some  point  and  continues  indefinitely  in  the  infinite  horizon  case.  Analytical 
issues  such  as  existence  and  convergence  of  the  infinite  horizon  cost,  and  optimality 
of  stationary  policies  have  explored  in  many  papers  such  as  Federgruen  and  Zipkin 
(1986a,  b),  DeCroix  and  Arreola-Risa  (1998),  and  Kapuscinski  and  Tayur  (1998). 
Especially,  Federgruen  and  Zipkin  (1987b)  show  that,  under  certain  assumptions, 
every  infinite  horizon  policy  has  finite  expected  discounted  cost  in  single-product 
case.  DeCroix  and  Arreola-Risa  (1998)  argue  that  every  finite  expected  discounted 
cost  in  multi-item  case  is  bounded  above,  and  show  the  existence  of  the  optimal  cost 
in  infinite  horizon  cases.  We  do  not  explore  detailed  analytical  issues  in  this  study. 
Instead,  we  show  the  convexity  of  the  expected  cost  function  for  the  infinite  horizon 
case,  and  determine  the  optimal  replenishment  policy. 

Let  is  the  discounted  factor,  where  0 < ft  < 1.  Then  the  finite  horizon  cost 
function  in  period  n is 

Gn(Sn,Qn)  = gn(Sn,Qn)  +PE  [G*+1(Sn+1)] , n=l,...,m  (2.19) 

where 


G“(S")  = n .,r  „ [S»(S».0»)+^[G;+1(S„+1)]] 

V»,n-  1 V^tit Qn m3 

and  G*m+l(Sm+l)  = 0,  1T  = (1, ...,  1). 

Define  Zn(S'n,Q'n)  = Gm_n+1(Sm_n+1,  Qm+n-i)  and  zn(S'n,Q'J  = 
gm-n+i{Sm-n+i,Qm+n-i)-  Then  Gn(Sn,Qn ) can  be  expressed  as  follows:  For  n = 


1,  ■■■>  m 
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Zn(S'n,  Q'n)  = zn(S'n,  Q'n)  + 0E  [Z^S'n  + Q'n-  Dn)]  (2.20) 

where  Dn  — ( Di>n , D^n)  is  a vector  of  random  demands  and  Zq(S'0 ) = 0. 

Let  m = oo  and  n — » oo.  Then  the  second  term  in  the  right  hand  side  in  (2.20)  is 
iterated  recursively  for  indefinite  number  of  periods  in  order  to  compute  the  infinite 
horizon  expected  cost. 

We  now  show  the  convexity  of  Zn(S'n,Q'n ) as  n -*  oo.  Following  the  arguments 
of  DeCroix  and  Arreola-Risa  (1998),  every  Z*(S'n)  is  bounded  above  and  there  exists 
a function  Z*(S')  with 

Z-(S')  = lim  Zn'(S'„) 

n-KX3 

since  Z*(S'n ) is  nonnegative  and  thus  is  increasing  in  n for  fixed  S'n.  That  is,  the 
infinite  horizon  cost  is  convergent  and  is  well  defined.  Analogous  to  above  notation, 
let  Z(S')  = lim^oo  Zn(S'n,  Q'n ) and  let  Q £°  be  the  optimal  solution  in  infinite  horizon 
case,  starting  in  period  n.  We  then  have  the  following  proposition: 

Proposition  2.7  The  expected  holding  and  penalty  cost  for  periods  n through  oo, 
Z(S')  is  convex  in  replenishment  sizes,  Q'n  where  Q'n  = (Q ij7l, ...,  QitU). 

Proof  : It  is  straightforward  to  derive  the  convexity  of  Z(S').  From  (2.20),  the 
optimal  finite  horizon  cost  is 


= Zn(S'„, <?„)  + PE  {ZU(S'»  +Q'„~  A.)]  , n = 1 m 
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and  when  m = oo  and  n -*  oo,  we  have 

Z-(S')  = z„(s;,Q”)  + PE  [z;_,(s;  + <?»  - d„)] 

As  we  argued  on  Theorem  2.1,  the  recursive  term  or  the  future  cost,  Z*_t,  is  always 
independent  of  decision  variables  whether  n is  a finite  or  goes  to  infinity.  Therefore, 
the  corresponding  term  is  vanished  in  the  first  and  second  partial  derivatives  with 
respect  to  Q'n  — Q £°,  and  our  task  reduces  to  derive  the  convexity  of  the  single-period 
cost,  zn(S'n,Q £°),  which  is  obvious  from  the  previous  result.  Therefore  Z*(S')  is 
convex  and  Z*(S')  ->  oo  as  ||S'||  — > oo.  As  a result,  Z(S')  is  convex  and  Z(S')  — > oo 
as  HS'll  — > oo  □ 

DeCroix  and  Arreola-Risa  (1998)  showed  the  convexity  of  the  infinite  horizon 
discounted  cost  in  multi-item,  periodic  review  inventory  system  under  the  assump- 
tion of  homogeneous  products.  They  also  argued  that,  in  case  of  non-homogeneous 
products,  they  could  not  provide  a practical  approach  for  solving  the  problem.  The 
result  of  Proposition  2.7  is  attractive  because,  in  this  study,  we  did  not  restrict  our 
assumption  to  the  homogeneous  product  case.  We  also  can  notice  that  our  method 
to  derive  the  convexity  of  the  expected  cost  function  in  either  a finite  or  an  infinite 
horizon  case  is  straightforward  and  easy  to  understand.  The  following  proposition 
provide  the  optimal  solution  of  the  infinite  horizon  case. 

Proposition  2.8  The  optimal  replenishment  policy  for  item  i in  the  infinite  horizon 


case  is 
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«°  = I-'*-—1.2-- 


where  p,  is  a Lagrange  Multiplier,  Fi<n  is  the  cumulative  distribution  function  of  a 
random  demand  DijTl,  and 


u = arq  min 

r y m>o 


V F~l 

/ -/  i,n 


L i=i 


i ( Pi~  V 
hi+Pi 


<V0 


Proof:  Define  L = Gn{h,n,  Qi,n, ...,  Qi,n)  + p (e!=i  ^,r»  + E!=i  Qi,n  - V&). 
Then,  from  (2.20),  we  get 


dL 


aG„+  S 


dQi,  n dQijU  dQi<n 

d 

-9n  + P 


mE  '<■" + E Q*  ~ v>) 


i=l  i=l 


dQi,n 

= {hi  + Pi) Fi^n{Ii^n  + Qi,n ) Pi  “t~  P — 0 


(2.21) 


Solving  (2.21)  and  the  KKT  condition:  p($2\=i  k,n  + Y!i=i  Qi,n  ~ Vo)  = 0,  we  get 


Q?n  = Fi 


i f Pi  + F 


oo 


v hi+Pi 
i 


- I, 


OO  


p = 


aromin 

M>0 


□ 


2.6  Conclusions 

In  this  chapter,  we  discuss  a stochastic  multi-item  inventory  model  for  the  case 
of  equal  replenishment  intervals  with  a warehouse-capacity  constraint.  We  com- 
pletely characterize  the  inventory  model  by  using  dynamic  programming.  We  use 
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induction  to  derive  the  convexity  of  objective  cost  function  and  show  that  a myopic 
replenishment  policy  is  optimal  for  both  a finite  and  an  infinite  planning  horizon. 
Furthermore,  we  notice  that  the  optimal  policy  in  each  period  is  easy  to  calculate 
due  to  the  independence  of  the  one  period  expected  cost  over  periods. 


CHAPTER  3 

UNEQUAL  REPLENISHMENT  INTERVALS 
3.1  Introduction 

In  this  chapter,  we  extend  our  study  to  a case  where  individual  items  have  differ- 
ent replenishment  schedules.  By  relaxing  the  assumption  of  identical  replenishment 
schedules  for  items  we  generalize  the  multi-item,  periodic-review  inventory  model 
described  in  the  previous  chapter.  Items  are  replenished  periodically  by  manufactur- 
ers, but  individual  items  have  different  replenishment  schedules.  Thus,  the  simple 
periodic-review  inventory  policy  is  not  valid  anymore  in  this  case.  This  assump- 
tion can  be  partly  explained  as  follows:  The  transportation  company  that  controls 
the  entire  distribution  system  can  efficiently  utilize  its  transportation  resources  and 
procure  its  transportation  services  at  lower  cost  by  assigning  a different  shipment 
schedule  to  each  product.  A warehouse  supplies  products  from  manufacturers  to  a 
population  of  customers.  We  first  state  the  model  assumptions  and  then  formulate 
the  expected  cost  function  to  determine  an  optimal  replenishment  policy 

3.2  Model  Assumptions  and  Formulation 
Customer  demands  for  items  are  random  variables,  and  backlogging  is  allowed. 
The  distribution  of  demand  for  each  product  is  assumed  to  be  stationary.  We  assume 
zero  replenishment  leadtime  between  the  warehouse  and  the  manufacturers,  and  all 
replenishment  orders  are  received  at  the  beginning  of  the  following  replenishment 
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period.  A finite  planning  horizon  is  considered.  We  use  the  following  notation  to 
describe  this  capacity  constrained  multi-item  inventory  model  with  different  replen- 
ishment schedules: 


Parameters 


l : number  of  items, 

m : total  number  of  time  instants  that  replenishments  occur, 

Tj  : time  instant  when  the  j th  replenishment  occurs,  j = 1, 

Tj  : length  of  replenishment  interval  for  item  i, 

T : length  of  the  time  horizon, 

V0  : fixed  storage  capacity  of  the  warehouse, 

hi  : unit  cost  of  holding  inventory  for  one  replenishment  period  for  item  i, 

Pi  : unit  penalty  cost  of  item  i for  demand  backordered  for  one  replenishment 
period, 


Demands 


di,t  '• 

Di>w 


demand  of  item  i during  time  (t,t  + 1),  a random  variable, 
where  t is  the  index  of  the  time  period,  t = 0, 1,  ...,T, 

: cumulative  demand  of  item  i over  w time  periods,  where  DiiW  = 
V t+w  d- 

/Ds=t+1 


State  and  decision  variables 


Iij  : inventory  level  of  item  i at  time  Tj  before  replenishment,  where  are 
given, 

Qij  : replenishment  quantity  of  item  i at  time  Tj. 


We  define  Rj  = {pij,  Pkjj}  as  the  set  of  items  replenished  at  time  Tj,  where 
kj  (<  l)  is  the  number  of  items  replenished  at  Tj  and  Pij  denotes  the  index  of  the  ith 
item.  We  also  define  Qj  = {QPlJ,  Qpkj,j)  as  a vector  of  decision  variables,  where 
Qij,  is  the  replenishment  quantity  of  item  i at  time  Tj,  where  i 6 Rj  . 

Since  we  allow  different  replenishment  schedule  for  each  item  (see  Figure  3.1),  it 
is  necessary  to  know  the  current  inventory  level  of  each  item  to  check  the  capacity 
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constraint  when  new  shipments  from  manufacturers  occur  at  the  warehouse.  Thus, 
determining  the  replenishment  size  of  an  item  depends  on  the  current  inventory  level 
of  all  other  items  as  well  as  the  warehouse  capacity.  We  notice  in  Figure  3.1  that 
the  replenishment  intervals  overlap  since  a replenishment  for  one  item  can  occur  in 
the  middle  of  the  replenishment  interval  for  other  items. 

We  first  define  Sj  = {Iij, ...,  Iij)  as  the  state  of  the  system  at  time  Tj  or  in 
stage  j.  For  a given  state  Sj,  the  current  inventory  level  of  the  system  in  stage  j is 
Ij  — Yl\=i  Ii,j>  which  can  be  expressed  recursively  as  follows: 


i 

Ij  = Ij- 1 + } ] QiJ—l  ~ y ] Di,Wj  (3-1) 

ieRj-i  i=l 

where  Wj  = Tj  — Tj_i.  We  also  determine  the  current  inventory  levels  of  individual 
items  recursively  using  the  state  in  the  previous  stage: 


IiJ 


lij— 1 T Qi,j— 1 Hi,Wji  if  i ^ Rj  — 1 

< 

Iij-i  — DitWj,  otherwise 


(3.2) 


We  define  gj(Sj,Qj ) as  the  expected  holding  and  penalty  costs  of  the  system 

in  stage  j,  given  initial  inventories  of  Sj  = (Rj,  and  replenishment  sizes  of 

Qj  = (QPlJ,  ...,QPk,  The  detailed  expression  of  gj(Sj,Qj)  will  be  shown  later. 

Denote  the  expected  holding  and  penalty  costs  for  stage  j through  m by  Gj(Sj,Qj ) 

and  let  Gj(Sj)  — min  Gj(Sj,Qj ) be  the  minimum  expected  costs  in  stage  j 
Qi,j  > i&Rj 

through  m,  given  the  system  starts  stage  j in  state  Sj  — (hj,  ■■■,  hj)  and  is  subject 
to  the  capacity  constraint,  Y?i= i I*J  + Q'j  < V0.  We  then  state  the  recursive 
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Figure  3.1:  A multi-item  inventory  system. 


equations  whose  solution  provides  an  optimal  policy  for  the  problem. 


Gj(Sj,Qi)  = Si(Sj.<3i)+£[Gj+1(Si+1)] 

roo  roc 

= 9j{Sj,Qj)+  ••■/  [G*j+l(Sj+i)]flJ(xi)---flJ(xl)dxl“-dxl 
Jo  Jo 

(3.3) 


where  is  the  probability  density  function  of  a random  demand,  Aj- 

We  eventually  want  to  find  the  optimal  cost  for  the  entire  planning  horizon, 
G*(/i),  which  can  be  calculated  recursively  by  solving  the  following  equation,  for 
j = 2, ...,  m, 
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Gj(Sj)  = 


min 


Qi,i:12i=i  Qij<Vo 


9j(Sj,Qj)  + f 

Jo 


poo 

/ [G*+1(5i+0]  ■ • • fij(xi)dxi  ■ • - dxi  (3.4) 

Jo 


We  now  express  the  single-period  expected  cost  in  any  stage  j as  follows: 
9i(Ss,Qj)  = ^ E + Qij  — Di}T.)+  + Pi{DitTj  — Iij  — Qi,j)+] 

ieRj 

roc 

— ^ + Qij  — Xi)+  -f-  Pi(xi  — /jj  — Qt,j)  ] fij(xi)dxi 

ieRj  Jo 

„ r 

= / y I / (Xi,j  "I-  Qij  ~ Xi)fij(Xi)dXi 
icB,  L Jxi=0 


+ 


poo 

p I ( xi  — Iij  — Qij)fij(xi)dxi 

J Xi  = Iij  +QiJ 


(3.5) 


where  fij  is  the  probability  density  function  of  the  demand  for  item  i in  period  j. 

Note  that  this  expected  cost  only  includes  those  items  that  are  replenished  at 
time  Tj.  The  reason  for  this  can  be  explained  as  follows:  for  simplicity,  we  only 
consider  the  first  two  items  in  Figure  3.1.  Since  item  2 is  scheduled  to  be  replenished 
at  time  Tj,  we  consider  one  replenishment  interval  of  the  item,  which  ranges  between 
Tj  and  Tj  + T2.  Then  the  total  cost  in  this  stage  is  the  sum  of  the  cost  for  item 
1 incurred  between  Tj  and  Tj  + T2  (say,  cost  1)  and  the  cost  for  item  2 incurred 
between  Tj  and  Tj  + T2  (say,  cost  2).  Cost  1,  however,  is  not  affected  by  the  current 
decision  variable,  Q2j,  whereas  the  corresponding  item  provides  information  about 
its  current  inventory  level  to  the  system  in  order  to  check  the  capacity  constraint, 
£<=i  Ji,Tj  + Q?  j < Vq.  That  is,  the  information  about  the  current  inventory  level  of 
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item  1 is  embedded  in  the  capacity  constraint.  How  about  the  cost  related  to  item 
1 (cost  1)?  As  we  mentioned  before,  cost  1 is  not  affected  by  the  decision  variable, 
Qi,j-  Thus,  when  we  determine  the  optimal  policy,  the  final  decision  is  not  affected 
by  whether  we  include  the  cost  in  the  objective  function  or  not. 

We  now  look  at  stage  (j  — 1),  which  runs  from  Tj_i  to  7j_i  +7\.  In  this  stage,  the 
decision  variable  is  Qij-i.  Note  that  the  cost  for  item  1 incurred  between  to 
Tj—i  + T\  includes  cost  1.  That  is,  cost  1 is  double  counted  in  periods  j and  ( j - 1). 
This  also  happens  to  item  2 (cost  2 is  counted  twice)  and  the  cost  is  not  affected  by 
Qij-i.  To  avoid  this  double  counting,  we  remove  one  from  where  it  is  not  necessary 
to  be  counted  (for  instance,  cost  1 in  stage  j and  cost  2 in  stage  ( j — 1)).  Once  we 
remove  it,  the  total  expected  holding  and  penalty  cost  in  stage  j only  includes  the 
cost  for  item  2 as  in  (3.5).  For  any  stage,  these  double-counted  costs  can  be  removed 
in  this  fashion,  and  finally,  the  expected  cost  for  the  entire  system  is  well  defined. 
For  more  than  two  items,  we  can  easily  extend  this  to  express  the  corresponding 
expected  cost  by  using  the  aforementioned  method. 

3.3  Solution  Approaches 

We  next  discuss  solution  approaches  to  the  dynamic  recursive  equation  (3.4), 
which  is  subject  to  the  capacity  constraint.  In  the  previous  chapter  we  established 
a solution  approach  for  minimizing  expected  cost  in  the  case  of  equal  replenishment 
intervals  (periodic- review  case).  We  also  found  that,  in  any  stage  j , the  optimal 
replenishment  quantity  was  independent  of  the  future  cost  incurred  in  that  stage. 
This  argument  is  generally  not  true  anymore  in  the  case  of  unequal  replenishment 
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intervals,  except  in  the  final  stage  of  the  planning  horizon  (stage  m in  our  case). 
To  show  this,  we  first  start  with  a simple  case  of  a two-item,  two-period  inventory 
system,  assuming  that  each  replenishment  occurs  in  time  periods  1 and  2. 

Two-item,  Two-period  Case 

Consider  the  warehouse  inventory  system  with  l = 2 items  and  m = 2 periods 
where  item  1 is  replenished  at  time  T\  = 1 and  item  2 is  replenished  at  time  r2  = 2. 
See  Figure  3.2.  We  use  same  notation  that  we  defined  earlier  in  this  chapter,  except 
for  defining  Qj  as  the  replenishment  quantity  in  stage  j,  for  simplicity. 

Given  that  J1( i and  /2)i  are  known,  the  inventory  level  in  period  (or  stage)  2 can 
be  expressed  as  follows: 


h, 2 — h, i + Qi  — rfi.i,  h,2  — h,i  ~ d2,i 


Note  that,  from  Figure  3.2,  the  expected  cost  in  period  2,  g2 (/^ 2,  /2)2,  Qi),  only 
includes  the  expected  cost  of  item  2 incurred  in  period  2,  whereas  the  expected  cost 
in  period  1,  <7i(/i,i,/2,i,Qi),  consists  of  the  expected  cost  of  item  1 incurred  in  one 
full  replenishment  period  from  time  1 to  T as  well  as  the  expected  cost  of  item  2 
incurred  in  period  1.  Then  the  expected  cost  for  period  2 through  the  end  of  planning 
horizon,  G2(ii>2,  /2>2,  Q2),  is  identical  to  </2(/ii2,  /2>2,  Q2)  and  can  be  stated  as  follows: 


G2(/ij2,  h,2,  Q2)  — <72(A,2,  T2,2j  Q2)  + E[G\(1 1)2  — di>2,  /2)2  + Qi  — d2>2)] 


— Ed2  2 [h2 (12,2  + Q2  — d2)2)+  +P2(d2)2  — /2)2  — <?2)+]  (3.6) 
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Item  1 


Item  2 
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| 

/1.1  V 

h,2  -i 

"" 
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J 

h,2  ~ 

1 

Figure  3.2:  A two-item,  two  period  inventory  system. 

(assuming  that  G\  = 0).  Applying  the  same  method  as  in  the  periodic-review  case, 
we  can  solve  (3.6)  to  determine  Q2,  which  becomes  a function  of  Q 1,  di,i  and  d2)i, 
as  follows: 


Q2(Qu  ^t,tj  ^2,1)  — Fr 


-1  f P2  P2{Qi  > ^1,1) 


= Fr 


2,2  ' h2+P2 

-1  f P2-  P*2{Qi,di,i) 


/4(Q  i)  di,i) 


2,2 


arq  mm 
/«>o 


= arq  mm 

H2>0 


h2  + p2  ) 
p-t  ( P2  P2 


2’2  Vh2  + 


— h,2 

~ {h,i  ~ d2>  1) 
+ h, 2 < Vo 


2 + P2 


F2  2 ( ) + /i,i  + Qi  — di  1 < V0 

’ V h2  + P2 ) 


or  simply 


Q2{Qi,di,i,d2,i)  = min 


P2 


h2  + p2 


^ » Vo  - (A,t  + Qt  - dt,i) 


72,i  + d2>i 


(3.7) 


where  Fjj  is  the  cdf  of  dy.  Note  that  p*  here  has  an  extra  term,  /1; 2,  compared  to 
that  in  the  periodic-review  case  in  the  previous  chapter,  (2.18).  /i>2  is  the  current 
inventory  level  of  item  1 which  is  not  replenished  in  period  2.  Since  all  items  are 
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replenished  at  the  same  time  in  the  periodic-review  case,  we  do  not  have  this  extra 
term  in  that  case.  Even  though  we  do  not  replenish  item  1 in  period  2,  we  still  need 
information  about  inventory  levels  of  this  item  to  check  the  capacity  constraint.  The 
expected  cost  function  for  period  1 through  the  end  of  planning  horizon  is  then 


9i(h,i,h,i,Qi)  + Edll  [G*2(h, l + Qi  - di,uh,i  - <£2,1)] 

ED,  [/ii(ii,i  + Qi  — D\)+  + p\(Di  — Iiti  — Qi)+] 

+-£^2,1  \hl{h,l  — d2,l)+  T P2 (^2,1  - ^2,l)  + ] 

+Edi,i  [Ed2,i  [Ed2,2  [£2  (-£2,1  + ^2,1  + Q2 (Qi)  *£1,1  > ^2,1)  - d2,2)+ 

+ Vi  (d2,i  — -£2,1  — ^2,1  — Q2(Qi,  diti,  o?2,i))+] ] ] (3-8) 


where  Di  = diti  + <£i)2.  The  first  term  in  the  right  hand  side  in  (3.8)  is  the  ex- 
pected cost  of  item  1 incurred  in  both  periods  1 and  2,  and  the  second  term  is  the 
expected  cost  of  item  2 incurred  in  period  1.  We  now  state  the  following  proposition. 


Proposition  3.1  The  expected  holding  and  penalty  cost  in  period  1 , Gi(/i,i,  h,i,  Qi), 
with  under  limited  capacity  is  a convex  function  of  the  replenishment  size,  Q\. 
Proof  : Let  L\  = Gi(/i,i,  Ii,i,  Qi)  + + -£2,1  + Qi  — Vo).  Taking  the  first  and 

second  partial  derivatives  on  Li  with  respect  to  Q\,  we  get  the  following: 

= dQi  -^.ij  Qi)  + (-£1,1  + h,2  + Qi  — Vo)] 

= \EDl  [hi(h,i  + Qi  — Dx)+  +pi(Di  — Iiti  — Qi)+]] 

Q 

+-^rEdl  1Edil  [Ed22  [h2  (£2,1  - d2,i  + Q2(Qi,  dip,  d2,i)  - d2,2)+ 

°Vi 

+P2  {d2,2  ~ h,i  + <£2,1  — Q2{Qi,  diti,  o?2,i))+] ] + V'X 
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— {hi  + pi)Fi(Iiti  + Qi)  — pi  + Hi 
d 


an  l f [^2  (A(Qi,diti)  - d2,2)+ 

C'Vl  J(h,l  Jd.2,2 

+P2  (c/2,2  — A(Qi,  dj.,l))  + ] /2,2/l,ldrf2,2drfl,l 


— (hi  +Pi)Fi(Iiti  + Qi)  — pi  + pi 


d 


dQ 


/i2  (A(Qi,  diti)  — G?2,2)  /2, 2*3^2, 2 

1 J d\ti  | J d2,2=0 
poo 

+ P2  (^2,2  “ ^(Ql,dl,l))  /2,2dd2,2 

•/  d2,2=A(Qi,di,i) 


/l.iddi,!  (3.9) 


(hi  + pi)-Fi(/i,i  + Qi)  — Pi  + Pi 

d ro° 

+Wi 


*F, 
' d 2,2 


r°°  rp2,2 

/ h2  (F2  2 — d2j2)  h,2dd2,2 

Jd\ti=a(Qi)  J d2,2=0 

r*OC 


poo 

+ / P2  (d2,2  - F^2  ) /2,s 

Jd2,2  = F2,l 


d d 


2UU2,2 


fi,idd 


1,1 


a ra(Qi)  rKQ  1) 

an~  / h2(b(Qi)  - d2t2)  f2j2dd2t2 

°Wl  J di,i=0  Jd2,2=0 


f 


+ / j?2  (o?2,2  — KQl))  f2,2dd2t2 

' d2,2=b{Qi) 


f i,iddi,i 


(3.10) 


where  P\  is  cd/of  D1;  /ii#  is  pdf  of  ditj,  A(Qi,dhi)  = min  , V0  - h,  2 = 


mm 


from  (3.7),  a(Qi)  = F^2  - V0  + JM  + Qx 
and  b(Qi ) = Vo  — A,i  — Qi  + di,i-  a(Qi)  and  b(Q  1)  in  (3.10)  are  obtained  as  follows: 
Prom  (3.9),  we  get 

Fd  (s&j)  . if  f;,1  (>£*)  <Vo-  hf 

Vq  — /i,2  , otherwise 


A(Qi,  dltl)  — < 


and  then 
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P2 


h2+P2 
P2 


and 


< Vo  — h, 2 

(fe+s)  £ y° “ (/u  + 0l  _ dw) 
di,i  > -F2T2  +P2)  ~ ^°  + ^1’1  + Q1  = °(<30 

fe(Qi)  = Vo  - h>2  = Vo  - /i,i  — Qi  + di,i 


Returning  back  to  (3.10),  we  have 


dLi 

<9Qi 


(/ii  + pi)Fi(/i)i  + Qi)  — pi  + pi 

»F, 

^2,2 


+ 


r^2,2 

/ ^2  (-^2,2  “ ^2,2)  /2,2dc?2,2 

^2,2=0 

poo 

+ I P2  (d'2,2  — ^2,2)  f 2, 2^2, 2 

Jd2,a=F^2 

rF2,i 

/ In  (J#  - ^2,2)  12,2^2,2 
J ^2.2=0 


/i,i(«(0i)) 


2^^2,2 


poo 

+ / P2  {d2,2  - F2,2  ) /2,: 

V d2(2=^*2,2 

ra{Qi)  a [ rb(Qi) 

+ / o^T-  / h2  (b(Qi)  — ^2,2)  /2,2dc?2,2 
Jdl,l=0  aWl  [Jd2, 2-0 

poo 

+ / P2  (^2,2  — b(Ql))  f2,2dd,2,2 

J d2,2=b(Qi) 

= (/ii  +pi)-P’i(A,i  + Qi)  — Pi  + Mi 


/l,l(o(Ql)) 


/i,iddi,i 


+ 


/■o(Qi)  ,0  [ /*6(Qi) 

/ ~ ^2>2)  ^2>: 
4dl,l=0  a(Vi  [4^2, 2=0 

poo 

+ / P2  (^2,2  — b(Qi))  /: 

^ d2.2=6(Ql  ) 


>dd' 


,2UU2,2 


2,2^2, 2 


/l.iddi,! 


' d2,2=b{Qi ) 

= (/ii +pi)i7’i(Zi,i  + Qi)  — Pi  + /^i 

ra(Qi)  rKQ\)  r°° 

+ / — h2f2,2&d2,2  + / P2f2,2&d2,2 

J di:i=Q  J d2,2=0  Jd2,2=b(Qi) 


/l,lddi,i 
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= (hi  + pi)Fi(ii,i  + Qi)  — Pi  + Pi 


r<Qi) 

rb(Q  1) 

rb(Qi) 

~h2  / 

h, 2^2,2  ~P2  /2,2dG?2,2  + P2 

J di,i  =0 

J ^2,2=0 

J d2,2=0 

/i,iddi,i 


= (hi  + pi)Fi(Iiti  + Qi)  — pi  + fj-i 

-a(Qi) 


d 2 

dQi2 


ra{W\) 

+P2'P1l,l(a(Ql))  ~ (^2  +P2)  / -^2,2 (^(Ql))^l,ldrfl(l 

Jdi,i=0 

(/ll  +Pl)/l,l(/l,l  + Ql) 


+ (h?  + P2) 


P2 


h2  + P2 . 

a(Qi) 


hMQi)) 


+ (/l2  + P2)  / /2,2(KQl))/l,l(a(Ql))^^l>l  — ^ ^ 

«/  di  1 =0 


(3.11) 


The  result  in  Proposition  3.1  shows  that  the  two-item,  two-period  inventory 
model  has  an  unique  optimal  replenishment  policy  and  the  solution  satisfying  the 
first-order  conditions  is  optimal.  This  solution  approach  would  work  in  rather  simple 
cases  with  small  number  of  items  or  periods.  When  the  number  of  items  or  peri- 
ods increase,  the  number  of  terms  to  be  evaluated  (e.g.  those  in  (3.10))  increases 
rapidly  and  it  would  be  impractical  to  solve  the  problem  using  the  optimal  solution 
approach.  This  leads  us  to  propose  heuristics  to  determine  replenishment  policies. 

3.4  Heuristic  Algorithms 

For  clarity  of  presentation,  we  first  consider  two-item  inventory  systems,  assuming 
that  only  one  item  is  replenished  in  a period.  Two  distinct  replenishment  schedules 
are  used  to  formulate  heuristics  and  we  then  extend  those  heuristics  to  more  general 


cases. 
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3.4.1  Heuristics  for  a Basic  Replenishment  Schedule 

The  first  replenishment  schedule  we  consider  is  where  item  1 is  replenished  in 
every  odd  period  and  item  2 is  replenished  in  every  even  period  until  the  end  of  the 
time  horizon  (see  Figure  3.3).  We  call  this  replenishment  schedule  Schedule  I. 

The  first  heuristic  (Heuristic  A)  is  motivated  from  Proposition  2.6  in  the  previous 
chapter  where  we  dealt  a problem  with  equal  replenishment  intervals.  Referring  to 
(2.17)  and  (2.18),  we  can  see  that  when  the  capacity  constraint  is  tight,  the  overall 
capacity  of  the  warehouse  in  any  period  j,  V0>  is  divided  into  two  terms  as  follows: 


( Pi  ~ Pi  A 
\h\  + pi  J 


( P2  ~ Pj\ 
\h2  +p2) 


= V0, 


Pj  > 0 


(3.12) 


Note  that  each  of  terms  on  the  left  hand  side  in  (3.12)  is  just  a function  of  unit 
holding  and  penalty  costs  and  the  demand  distribution  of  items  1 and  2,  respectively, 
and  those  values  are  given.  The  first  term  is  related  to  item  1 and  the  second  term 
is  related  to  item  2.  When  the  capacity  constraint  is  not  tight,  we  have 


F-i  (K-K) 

1j  Ui+pJ 


(P2  ~ Pj\ 

\ h-2  + P2  J 


<V0, 


Pj  0 


(3.13) 


We  now  apply  these  results  to  the  case  with  unequal  replenishment  intervals.  We 
assume  that  the  demand  of  each  item  is  identically  distributed  over  periods,  and 
thus  we  can  drop  subscript  j.  We  project  the  replenishment  schedule  of  item  1 and 
item  2 in  Figure  3.3  so  that  both  items  have  the  same  replenishment  schedule  (It  is 
possible  only  in  this  particular  case  because  the  length  of  replenishment  intervals  for 
each  item  is  identical.  We  consider  more  general  cases  later.)  We  then  consider  the 
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Figure  3.3:  An  inventory  system  with  unequal  replenishment  intervals  (Schedule  I). 

first  terms  in  both  (3.12)  and  (3.13)  as  the  capacity  of  item  1,  Vi,  and  the  second 
terms  as  the  capacity  of  item  2,  V2,  over  the  entire  planning  horizon.  That  is,  we 
split  the  overall  capacity,  Vo,  into  two  individual  capacities.  We  then  determine  Vi 
and  V2  by  solving  the  following: 


Find  n*  = min{/i  > 0} 


s.t.  Vi  + V2  < Vo 


Vi  = Ff 1 
V2  = Ff1 


P i-V 
hi  +Pi 


( P2~  P\ 

\h2  + p2J 


(3.14) 


Note  that  Fj  is  the  cumulative  distribution  function  of  item  i for  two  periods  in 
this  case.  Once  we  determine  those  individual  capacities  for  items,  we  return  back 
to  our  original  problem  and  assume  that  each  item  has  a separate  capacity  (see  Fig- 
ure 3.4)  The  replenishment  quantity  of  the  item  replenished  in  period  j is  limited  by 
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Item  1 


Item  2 


Figure  3.4:  An  inventory  policy  of  Heuristic  A. 

its  own  capacity.  That  is,  when  we  determine  the  replenishment  quantity  of  an  item 
in  period  j,  we  only  need  to  consider  the  corresponding  capacity  and  current  inven- 
tory level  of  that  item.  Once  each  item  has  its  own  capacity,  we  have  two  separate 
periodic-review  inventory  systems  and  thus  we  can  determine  the  respective  optimal 
replenishment  quantity  (for  the  separated  problem)  of  each  item  separately.  Note 
that  individual  capacities  depend  on  demands  and  unit  costs  and  thus  reflect  this 
information  when  determining  capacities.  This  heuristic  always  replenishes  item  i up 
to  V^,  i—  1, 2.  Referring  to  (3.13),  we  can  see  that  the  size  of  Vi  depends  on  whether 
the  overall  capacity  constraint  is  tight  (fi  > 0)  or  not  (//  = 0).  When  demand  rates 
are  not  high,  compared  to  the  warehouse  capacity,  the  overall  capacity  constraint 
is  not  tight  and  the  sum  of  two  individual  capacities  should  be  less  than  the  total 
capacity,  V0.  Then  the  heuristic  has  a replenishment  policy  which  replenishes  item  i 
up  to  the  unconstrained  optimal  inventory  level  or  V*  = Ff1  f J , assuming  that 
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item  i is  replenished  in  period  j,  and  the  policy  becomes  optimal  for  the  constraint 
case  as  well.  When  the  constraint  is  tight,  Heuristic  A replenishes  item  i up  to 
Vi  = Ft~l  where  Vi  + V2  = V0.  One  possible  disadvantage  of  this  heuristic 

is  that  the  warehouse  might  have  unused  capacity  for  one  item  while  the  other  item 
needs  a greater  replenishment  quantity  than  its  own  capacity  to  meet  demand.  In 
most  cases,  however,  this  situation  would  be  resolved  by  the  dependency  of  capacity 
on  demands  because  when  one  item  has  a higher  demand  the  corresponding  capacity 
would  be  increased  as  well.  Therefore,  overall,  this  capacity-separate  scheme  would 
be  plausible  even  though,  in  some  individual  periods,  the  aforementioned  disadvan- 
tages might  occur.  Finally,  each  item  is  replenished  up  to  its  own  capacity. 

The  second  heuristic  (Heuristic  B)  is  similar  to  the  first  one  in  the  sense  that  we 
use  a separate  capacity  for  each  item.  Instead  of  using  fixed  capacities,  Vi  and  V 2,  this 
heuristic  adjusts  those  capacities  according  to  the  availability  of  any  unused  capacity. 
This  adjustment  is  based  on  the  following:  Assuming  that  item  1 is  replenished  in 
period  j and  item  2 is  not,  the  inventory  level  of  item  2 is  less  than  its  own  capacity, 
V2,  because  the  inventory  level  has  been  decreasing  ever  since  the  last  replenishment 
of  item  2 (see  Figure  3.5).  Thus,  there  is  an  unused  space  for  item  2 at  the  beginning 
of  period  j and  we  want  to  allocate  at  least  a part  of  this  unused  capacity  to  item  1 so 
that  item  1 has  more  available  space  to  accommodate  its  demand.  Then  a question 
arises  immediatelyhow  much  we  can  increase  the  capacity  of  item  1.  Basically,  we 
increase  it  as  much  as  the  average  demand  of  item  1 incurred  between  time  j and 
(j  + 1)  so  that  the  capacity  of  item  1 eventually  is  reduced  back  to  Vf  at  the  beginning 
of  period  ( j + 1),  and  the  capacity  of  item  2 then  decreases  as  much  as  its  average 
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Figure  3.5:  An  inventory  policy  of  Heuristic  B. 

demand  for  one  period  as  long  as  the  overall  capacity  in  period  j does  not  exceed 
V0.  As  a result,  the  newly  adjusted  capacity  can  then  be  determined  by  solving  the 
following  equations 

V1+aE[Dl]  + V2-(l-a)E[D2]  = l/0 

aE[Dx]  - (1  - a)E[D2]  = 0 (3.15) 

where  a is  a adjustment  factor  and  solving  (3.15)  results 

Vx  = V\  + aE[Dx\ 

V'  = V2  - (1  - a)E[D2] 

E[D2 } 

a E[DX ) + E[D2] 


If  item  2 is  replenished  in  period  j,  we  get 
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V\  — (1  — a)E[D\]  + V2  + otE\D-i\  — Vo 


and 

^ = Vi  - (1  - a)E[Dx] 

V2  = V2  + 01.Ej.D2] 

g[Al 

£■[£>!]  + E[Lh] 

We  apply  this  scheme  in  any  period  by  increasing  the  capacity  of  the  item  re- 
plenished in  that  period  and  by  decreasing  that  of  the  other.  Once  the  capacity  is 
adjusted,  an  individual  item  is  replenished  up  to  its  inventory  level,  which  is  the 
minimum  of  its  adjusted  capacity  (Vfc),  the  unconstrained  optimal  inventory  level 
after  replenishment  (i^_1  )),  or  the  total  available  capacity  (V0  - hj),  assum- 

ing that  item  i is  planned  to  be  replenished  in  period  j , where  k = 1,2,  A;  / i.  The 
adjustment  of  capacity  wouldn’t  hurt  the  overall  performance  of  Heuristic  B because 
it  only  allocates  any  available  space  to  where  it  needs  and  does  not  force  item  2 
to  reduce  its  physical  inventory  level  due  to  the  capacity  reduction.  Even  though 
item  1 has  a capacity  increment,  it  adjusts  its  inventory  level  only  to  a direction  to 
improve  the  overall  performance.  Therefore,  this  heuristic  would  improve  the  overall 
performance  from  Heuristic  A. 

In  the  last  heuristic  (Heuristic  C),  we  consider  the  total  available  capacity,  V0  — 
I2j,  in  the  system  at  the  time  instant  of  replenishment  instead  of  using  separate 
capacity  for  each  item.  Referring  to  Figure  3.6,  item  i is  replenished  up  to  its  total 
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e>J 

Item  1 

'M 

Item  2 

/ / +1 


Figure  3.6:  An  inventory  policy  of  Heuristic  C. 

available  capacity,  or  the  unconstrained  optimal  inventory  level  after  replenishment, 
F~l  ^ /iH-pT ) > whichever  is  smaller.  Thus,  heuristic  C always  maintains  the  greatest 
capacity  for  replenishment  among  three.  Note  that  when  the  capacity  constraint  is 
not  tight,  all  three  heuristics  have  the  same  replenishment  policy  which  replenishes 
up  to  F~l  We  summarize  these  three  heuristics  in  the  following. 


Heuristic  A 

• Compute  the  separate  capacity  of  each  item  by  solving  the  following: 


Find  n*  — min{/i  > 0} 


Vi  + y2  < vQ 


Vi  = Ff1 
V2  = F2~l 


Pi  — M 
hi+Pi 


( P2~  P\ 
\h2+  p2) 


s.t. 
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• Determine  the  replenishment  quantity 


Qi,j  — hi  Ii,j  > ® 


Heuristic  B 

• Compute  the  separate  capacity  of  each  item  as  for  Heuristic  A. 

• Adjust  the  capacities  of  individual  items  based  on  their  demands  as  follows: 

• If  item  1 is  replenished  in  period  j, 

V/  = Vx+aE[Dx] 

V'2  = V2-{1  -a)E[D2] 

E[D2] 

a E[Di]  + E[D2] 

■ If  item  2 is  replenished  in  period  j, 


Vi  = Vi  — (1  — a)E[Di] 

V'  = V2  + aE[D2] 

E[DX] 

“ E[Di]  + E[D2] 

• Determine  the  replenishment  quantity 


Qu  - min 


Pi 


hi  + pi 


Vi,  V0 -I, 


k,j 


Iitj  , i,  k — 1)  2,  k ^ i 


Heuristic  C 

• Assuming  that  item  i is  replenished  in  period  j,  determine  the  available  capac- 
ity; Vq  — h,j,  i,  k = 1, 2,  k ^ i for  item  i. 
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• Determine  the  replenishment  quantity 


Qu  = min 


f: 


Pi 

hi  + pi 


Vo  - h 


> 


i,k  = 1,2 ,k  ^ i 


3.4.2  Heuristics  for  An  Alternative  Replenishment  Schedule 

We  now  consider  a more  general  replenishment  schedule  (Schedule  II).  Let  item  1 
be  replenished  in  periods  1,  3,  5,  7,  ...  and  item  2 be  replenished  in  periods  4,  8,  12, 
16,  ...  See  Figure  3.7.  To  determine  individual  capacities  of  items,  we  first  project 
the  replenishment  schedule  of  item  2 on  that  of  item  1 as  in  the  previous  case.  As 
we  can  notice,  the  length  of  replenishment  interval  of  each  item  is  not  equal.  We 
then  find  the  least  common  multiple  of  individual  interval  lengths,  which  is  four  in 
this  case  (we  call  this  the  length  of  one  replenishment  cycle).  Since  item  1 has  two 
replenishments  in  one  replenishment  cycle  while  item  2 has  one,  the  expected  cost 
for  one  replenishment  cycle  can  be  rewritten  as  follows: 


E[2hi(yi  — di)+  + 2p\(d\  — yi)+  + /i2  (2/2  — d2)+  + p2(d2  — 2/2)+] 


where  j/j  is  replenish  up  to  level  and  is  equivalent  to  the  inventory  level  of  item  i 
after  replenishment.  By  minimizing  the  expected  cost  with  respect  to  under  the 
warehouse  constraint,  the  following  is  obvious  from  the  previous  derivation 


V*i 

V*2 

P* 


F-i  ( 2?i  ~ P*  \ 

1 \2(hl+Pl)J’ 
F-  ( jZ-£) 

\h2Fp2J 


arq  min 
n>o 


f 2pi  — p \ 
\2{hi  +Pi)J 


+ FP  (tt—E) 

\h2  +P2J 
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Figure  3.7:  An  inventory  system  with  unequal  replenishment  intervals  (Schedule  II). 

Using  the  same  approach  as  in  the  previous  section,  we  determine  individual  capacity 
of  each  item  by  solving  the  following 


Find  n*  = min{/x  > 0} 

s.t.  Vi+V2<  Vo 


Vi  = Ff1 
V2  = FT1 


( - n \ 

\2(hi  +Pi)J 

( P2~P\ 

\ /i2  + p2  / 


and  Heuristic  A has  a replenishment  policy  which  replenishes  item  i up  to  Vj. 

Heuristic  C remains  the  same  as  in  the  previous  case.  In  other  words,  the  heuristic 
replenishes  item  i up  to  its  total  available  capacity,  or  the  unconstrained  optimal 
inventory  level  after  replenishment,  Ff1  (^hp_lp. ) , whichever  is  smaller. 

For  Heuristic  B,  we  first  take  a look  at  Figure  3.7.  In  periods  j =1,  5,  9,...,  we  can 
see  that  no  replenishment  occurs  in  the  following  periods  (periods  2,6,10,...).  Hence, 
we  increase  the  capacity  of  item  1 as  much  as  the  average  demand  of  item  1 in  the 
current  and  following  periods  (periods  j and  j + 1),  so  that  the  adjusted  capacity 
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is  eventually  reduced  back  to  its  original  capacity  after  two  periods.  For  item  2,  we 
decrease  its  capacity  as  much  as  the  average  demand  between  its  last  replenishment 
period  (periods  0,  4,  8,  ...)  and  the  current  period,  which  is  equal  to  E[D2\.  We  then 
can  determine  the  adjusted  capacity  of  each  item  as  follows: 

Vi  + aE[2Di\  + V2  - (1  - a)E[D2\  = V0 

aE[2Dl]-(l-a)E[D2]  = 0 (3.16) 

where  a is  a adjustment  factor  and  solving  (3.16)  results 

V{  = V1  + 2aE[Dl ] 

V'  = V2-(l-a)E[D2] 

E[D2] 

a 2E[Dl]  + E[D2] 

In  periods  3,  7,  11,...,  a replenishment  occurs  immediately  in  the  following  periods 
4,  8,  12,...  and  item  2 was  not  replenished  for  the  three  preceding  periods.  Therefore, 
we  increase  the  capacity  of  item  1 as  much  as  E[D\]  and  decrease  the  capacity  of 
item  2 by  3E[D2]. 

Vi  + a£[A]  + V2  - (1  - a)E[3D2 } = V0 
ctE[Dx)  - (1  - a)E[3D2]  = 0 
=►  Vi  =Vi  + aE[Di] 

y-=v2-(i  a - 
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In  periods  4,  8,  12,...,  replenishment  occurs  immediately  in  the  following  periods 
5,  9,  13,...  while  item  1 is  not  replenished  for  only  one  preceding  period.  We  thus 
increase  the  capacity  of  item  2 by  E[D2]  and  decrease  the  capacity  of  item  1 up  to 
E[D2}. 

Vi  — (1  — a)E[Di]  + V2  + aE[D2]  = Vo 
— (1  — a)E[Di]  -F  cxE^D^  — 0 


V;  = Vi  - (1  - o^ElD,] 

V2  = V2  + aE[D2] 

E[D  i] 

“ E[DX}  + E[D2) 

To  combine  these  three  cases,  let  aj  be  the  number  of  periods  between  the  jth 
and  ( j + l)th  replenishment  in  the  system.  For  example,  in  Schedule  II,  oi  = 2,  a3  = 
1,04  = 1)  05  = 2,  so  on.  Assuming  that  item  i is  not  replenished  in  period  j,  we 
define  eitj  as  the  number  of  preceding  periods  since  the  last  replenishment  of  item  i. 
For  instance,  e2,i  = 1,  e3)2  = 2,  64,1  = 1 so  on.  Then  we  have  the  following: 

• If  item  1 is  replenished  in  period  j, 

V/  = Vi  + a E[aj  Di] 

V2  = V2  - (1  - a)E[e2>j  D2] 

= e2j-  E[D2\ 

aj  + e2j  E[D2\ 
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• If  item  2 is  replenished  in  period  j, 


v;  = V,  - (1  -o)£[«uD,] 


V'j  = Vi  + 02j,  a = 


eij  E[Di) 


ei ,j  E[D i]  + aj  E[Di ] 


Once  the  adjusted  capacities  are  determined,  Heuristic  B replenishes  item  i up 
to  its  capacity  (V*),  the  unconstrained  optimal  inventory  level  after  replenishment 
w1  (sSs)).  or  total  available  capacity  (Vo  — hj,  k — 1, 2,  k / i)  whichever  is 
smaller.  We  notice  that  this  is  a general  version  of  Heuristic  B for  any  two-item  case, 
assuming  that  no  simultaneous  replenishment  occurs  in  a period. 

3.4.3  Heuristics  for  Multi-Item  Cases 

Based  on  the  results  from  the  two-item  cases,  we  can  extend  these  Heuristics 
to  multi-  item  cases,  assuming  that  there  is  no  simultaneous  replenishment  in  one 
period.  For  Heuristic  A,  let  l be  the  number  of  items  in  the  system  and  b be  the  least 
common  multiple  of  individual  interval  length  (or  the  length  of  one  replenishment 
cycle).  If  Tj  represents  the  replenishment  interval  length  of  item  i as  defined  earlier 
in  this  chapter,  6,  = &/T)  represents  the  number  of  replenishment  of  item  i in  one 
cycle.  Then  the  expected  cost  for  one  replenishment  cycle  is 


y ^E  [ft*  (hi(yi  - di)+  +pi(di  - ?/i)+)] 


and  the  following  is  immediate 
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/x  = arg  mm 
m>o 


E*rl 


i=  1 


biPi  ~ li 

M(hi+Pi ) 


< Vo 


The  capacity  of  each  item  can  then  be  obtained  by  solving  the  following: 


Find 

s.t. 


p*  = min{/x  > 0} 


J2Vi<VQ 


i= 1 


Vi  = F[ 


i ( kpi  - p 
bi{hi  +pi) 


, i = 1,2, I 


and  Heuristic  A replenishes  item  i up  to  Vj. 

It  is  straightforward  to  extend  Heuristic  B to  multi-item  cases  from  two-item 
cases.  Without  loss  of  generality,  assume  that  there  are  l > 2 and  item  1 is  replen- 
ished in  period  j.  Then  we  can  determine  the  adjusted  capacities  of  items  by  solving 
the  following  equations: 

i 

V\  + aE[aj  A]  + (V^  — (1  — a)E[eij  Di] ) = V0 

i=2 

/ 

aE[aj  Dx]  - - a)E[eid  A]  = 0 

i= 2 

and  we  get 


Vi  — Vi  T ocE]a,j  D\] 

V'{  = Vi  - (1  - a)E[ei>j  A],  * = 2,3,...,Z 

a _ Yli=2ei,j  F[Pj] 

aj  E[Di]  +X^!=2ei,j  ^[A] 
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The  replenishment  policy  for  Heuristic  B remains  the  same  as  in  the  two-item 
case.  No  change  is  made  for  Heuristic  C.  We  now  summarize  these  three  heuristics 
as  follows: 


Heuristic  A 

• Compute  the  separate  capacity  of  each  item  by  solving  the  following  problem: 


Find 

s.t. 


/i*  = min{/i  > 0} 


£><<  Vo 


i= 1 


Vi  = Ft 


-l 


biPi-y 

M(hi  + Pi ) 


, i — 1,2, I 


• Determine  the  replenishment  quantity 


Qi.j  — V ! * — 1>2,  ■■■l 


Heuristic  B 

• Compute  the  separate  capacity  of  each  item  as  for  Heuristic  A. 

• Determine  the  number  of  periods  between  jth  and  (j  + l)th  replenishment  in 
the  system,  aj. 

• Assuming  that  item  1 is  replenished  in  period  j,  determine  the  number  of 
preceding  periods  of  item  i since  its  last  replenishment,  eij,  i = 2, 3, ...,  I . 

• Adjust  the  capacities  of  individual  items  based  on  their  demands  as  follows: 


— Vi  T oiE\flj  Di] 
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V[  = Vi  - (1  - a)E[eid  Di]  , z = 2,3,..i 

„ _ Ei=2  ei  J E[Dj\ 


aj  ElDi}  + EL  ehi  ElDi\ 


• Determine  the  replenishment  quantity 


Q i,j  = min  Fl  1 


Heuristic  C 


• Assuming  that  item  1 is  replenished  in  period  j , determine  the  available  ca- 


3.4.4  Heuristics  for  Simultaneous  Replenishment  Case 

The  heuristics  that  we  have  developed  so  far  are  limited  for  the  cases  that  only 
a single  item  is  replenished  in  one  period.  In  real  warehouse  systems,  however,  it  is 
more  likely  that  multiple  items  are  replenished  at  the  same  time.  The  main  issue  in 
this  case  is  how  to  allocate  the  available  capacity  to  those  items  replenished  simul- 
taneously. For  Heuristic  A,  the  process  to  determine  individual  capacities  remains 
the  same  as  the  single  item  replenishment  cases  since  the  heuristic  depends  on  the 
number  of  items  in  the  system,  but  not  on  the  number  of  items  replenished  simul- 


pacity, Vo  - E!=2^j>  for  item  1- 


• Determine  the  replenishment  quantity 


taneously.  This  is  true  because,  in  Heuristic  A,  we  determine  individual  capacity  of 
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each  item  under  the  assumption  that  all  items  are  replenished  at  the  same  time  and 
thus  it  does  not  affected  by  whether  an  item  is  replenished  or  not  in  a specific  time 
instant. 

To  formulate  the  second  heuristic,  let  J = {1,2, ...,/}  be  the  index  set  of  items 
in  the  system.  We  also  define  Rj  C J as  the  index  set  of  items  replenished  in  period 
j.  Similar  to  the  single  replenishment  case,  we  can  determine  the  adjusted  capacity 
of  each  item  by  solving  the  following  equations: 

^2(Vi  + aE[ajDi])+  ]T  (Vk  - (1  - a)E[ekJ  Dk})  = V0 

ieRj  keJ-Rj 

aE[aj  Di]  - ^ (i  - a)E[ekJ  Dk } = 0 

i€Rj  k£j-Rj 

and  we  get 

V-  ~ Vi  + aE[(ij  Di],  i G Rj 

Vk  =Vk  — (1  - a)E[ek}j  Dk],  k G J - Rj 

_ Y^kej-Rj  ek,j  E[Dk] 

YlieRj  ai  E[Di\  + YlkeJ-Rj  E[Dk\ 

We  observed,  in  the  single-item  replenishment  case,  that  Heuristic  C had  a re- 
plenishment policy  which  replenishes  a single  item  up  to  the  total  available  capacity 
in  a period.  We  follow  the  same  replenishment  policy  in  this  case.  However,  be- 
cause of  simultaneous  replenishment,  we  need  to  allocate  the  total  available  capacity, 
Vo  - YlkeJ-Rj  individual  items  that  are  replenished  in  the  same  period.  We 

first  determine  individual  capacities  of  these  items,  as  for  the  case  of  Heuristic  A. 
Note  that,  in  this  case,  we  only  compute  the  capacities  of  items  replenished  in  a 
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period  (say,  period  j)  instead  of  determining  the  capacities  of  all  items.  The  related 
expected  cost  in  period  j is 


^2  E [bi  (hi{yi  - di)+  +pi(di  - yi)+)] 

i£Rj 

where  bi  is  the  number  of  replenishments  of  item  i in  one  replenishment  cycle.  The 
optimal  order-up-to  level  of  item  i,i  6 Rj  under  the  warehouse  constraint  is  then 


y\  = Fi 


-i 


biPi  ~ p* 
bi{hi+Pi) 


p — arq  mm 
^ ^>o 


ieRj 


, i (E  Rj, 

( kpi  - p 


t)+  E 

l> ' keJ-Rj 


Mk+Pi ). 

and  the  capacity  of  each  item  can  be  obtained  by  solving  the  following  equations: 


Find  p*  = min{/i  > 0} 


s.t. 


Ev'<+  E 


i£Rj 

Vi  = Frl 


x £J—Rj 

( bjPi-p  \ 
\bi(hi  +Pi)J  ’ 


i e Rj 


Once  the  individual  capacities  are  computed  and  assigned  to  each  replenished  item, 
we  look  into  any  unassigned  storage.  Since  YlieRj  — Vq,  there  might  be  a certain 
amount  of  unassigned  space  in  the  warehouse  (we  call  thisleft-over  capacity).  For 
Heuristic  C,  we  use  up  all  available  capacity  for  replenishment.  Thus,  when  we 
replenish  items,  we  use  this  left-over  capacity  as  well  as  their  own  capacity.  For 
example,  when  we  replenish  the  first  item  (say,  item  ii),  the  current  left-over  capacity 
is  L\ o = Vo  — YlieRj  Vi~Yli'eJ-Rj  ,j  and  the  item  is  replenished  up  to  the  minimum 
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of  LV0  + Vjj  and  F j"1  yh  p'^p  j.  In  case  that  F^1  is  smaller,  the  left-over  capacity 
for  the  next  item  (say,  item  i2)  is  LVi  = LV0  + Vil  — F^1.  Otherwise,  there  is 
no  left-over  capacity  available  since  the  item  is  replenished  up  to  LV0  + V*,.  In 
other  words,  the  left-over  capacity  for  item  i2  is  LV2  — max{0,  LV0  + V*,  — F^1}, 
and  the  item  is  replenished  up  to  min{LV2  -I-  Vi2,  F^1}.  We  continue  this  process 
until  all  items  are  replenished.  Note  that  when  only  one  item  is  replenished  in 
period  j,  then  LV0  + Vix  = V0  - Yli'eJ-Rj  h'  ,j  an<^  the  item  replenished  up  to 
min{V0  - Ysi'ej-Rj  h' F^1},  which  is  identical  to  the  single-item  replenishment 
case  established  in  the  previous  section.  We  summarize  these  three  heuristics  in  the 
following: 


Heuristic  A 

• Compute  the  separate  capacity  of  each  item  by  solving  the  following  equations: 


Find 

s.t. 


H*  = min{/^  > 0} 


i=l 


Vi 


( hpi  - n \ 

\bi(hi  +Pi)J  ’ 


i — 1,2, ...,  I 


• Determine  the  replenishment  quantity 


(3.17) 


Qi,j  hi  Ii,j  ) * ^ Fj 

Heuristic  B 

• Compute  the  separate  capacity  of  each  item,  Vi,  as  for  Heuristic  A. 

• Adjust  the  capacities  of  individual  items  based  on  their  demands  as  follows: 
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• If  item  i is  replenished  in  period  j, 

V-  = Vi  + aE[aj  Di],  i G Rj 

where 

_ YlkeJ-Rj  ek,j  E[Dk] 

SieK,  aj  E[Di\  + YlkeJ-Rj  ek,j  E[Dk] 

• If  item  i is  not  replenished  in  period  j, 


V-  =Vi-{l-  a)E[eid  A],  ieJ-  Rj 


• Determine  the  replenishment  quantity 


Qi,j  = min 


l J 


k£j-Rj 


-l 


) 


i e Rj 


Heuristic  C 

Let  Rj  = ~;ik}  where  k is  the  number  of  items  replenished  in  period  j 

and  ik  is  the  index  of  the  kth  item  replenished  in  the  period. 

1.  Determine  the  separate  capacity  of  each  replenished  item  by  solving  the  fol- 
lowing equations: 


Find  p*  = min{^  > 0} 


s.t. 


I>+  £ <V4 


ieRj 

Vi  = Fr 


i'eJ-Rj 

I ( biPi  - p \ 
\bi{hi  + pi)J  ’ 


i € A 


2.  Compute  any  left-over  capacity  after  the  capacity  allocation  to  each  replenished 


item. 


70 


m , = v0-J>-  Y, 

i£Rj  i'zJ-Rj 

3.  Replenish  item  ii  up  to 

min{LVo  + Vii,  Fr1} 

4.  Compute  any  left-over  capacity  after  replenishing  item  ii. 

LV i = max{0,  LV0  + - F~r} 

5.  Repeat  3 & 4 for  item  i2,  ...,ik 

3.5  Numerical  Results 

This  section  presents  numerical  results  illustrating  the  performance  of  the  heuris- 
tics described  in  the  previous  section.  Three  different  sets  of  replenishment  schedules 
are  considered.  The  first  two  sets  are  Schedules  I and  II  which  we  mentioned  before. 
The  third  set  of  schedules  (Schedule  III)  is  for  the  simultaneous  replenishment  case. 
For  Schedule  III,  item  1 is  replenished  in  periods  1,  3,  5,  7,  ...  , item  2 in  periods  1, 
4,  7,  11,  ...  and  item  3 in  periods  2,  5,  8,  11,  ...  We  run  a number  of  sample  problems 
for  each  schedule  by  varying  input  parameters.  Three  different  levels  of  unit  cost 
parameters  are  chosen  to  be  hi,pi  = 2,5,  i = 1,2,3.  Demand  for  each  item  follows 
an  exponential  distribution  with  parameter  A i and  we  use  three  different  levels  of 
demand  parameters  for  each  item  (Aj=10,  30,  50).  Since  each  item  for  Schedule  I is 
replenished  in  every  other  period,  the  actual  demand  for  one  replenishment  period 
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has  an  Erlang  distribution  with  parameters  fi  — 2 and  A;,  i = 1,2  and  the  cumulative 
distribution  function  of  the  demand  is  then  as  follows: 


For  Schedule  II,  item  2 is  replenished  in  every  fourth  period  and  thus  its  demand 
has  an  Erlang  distribution  with  parameters  /3  = 4 and  A2  whose  cdf  is 


and  for  Schedule  III,  items  2 & 3 are  replenished  in  every  third  period  so  that  the 
respective  demand  for  both  items  have  an  Erlang  distribution  with  parameters  /3  = 3 
and  A*,  i = 2, 3 whose  cdf  is 


Unfortunately,  in  many  cases,  the  inverse  of  is  intractable,  so  we  employ  a numer- 
ical approximation  to  compute  Ft~l. 

The  initial  inventory  level  of  each  item  is  set  to  be  /g  1,  /2,i,  h,i  = 10,  respectively. 
The  choice  of  initial  inventory  levels  should  not  have  a major  impact  on  the  decision  of 
replenishment  policy  in  the  long  run.  The  system  is  simulated  for  m = 50  periods  for 
each  of  the  three  heuristic  policies  using  the  same  demand  stream  and  the  respective 
overall  expected  cost  is  recorded.  These  50-period  simulations  are  repeated  5,000 


times. 
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It  would  be  ideal  to  compare  the  expected  cost  of  each  heuristic  policy  to  the 
expected  cost  of  the  optimal  policy.  Unfortunately,  it  is  impractical  to  obtain  dy- 
namic programming  solutions  except  for  small  sized  problems  with  discrete  demands 
and  thus  such  a comparison  would  not  be  possible.  Even  though  we  do  not  try 
drawing  any  conclusions  from  the  results  obtained  by  running  small  sized  problems, 
we  might  get  insight  from  the  results  to  provide  some  measure  of  the  effectiveness 
of  the  heuristics  developed  in  this  study  to  the  optimal  policy.  For  this  purpose,  we 
run  a two-item,  two-period  problem.  Note  that  we  have  already  shown,  in  the  previ- 
ous section,  that  the  expected  cost  of  this  problem  was  convex  in  the  replenishment 
quantity.  Furthermore,  referring  to  (3.11),  the  optimal  solution  can  be  obtained  by 
solving  the  following  first-order  conditions: 


(hi  + pi)Fi(Iiti  + Q i)  — Pi  + pi 


^i(h,i  + h,i  + Qi  — Vq)  — 0 


where 


F2j2(x)  = 1 - e 


b(Q  i)  — Vq  — h,i  — Qi  + diti 
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We  solve  these  equations  using  numerical  approximation  since  the  solution  is  not 
tractable.  Tables  3.1  ~ 3.4  contain  the  average  costs  (over  the  planning  horizon)  of 
each  heuristic  policy  as  well  as  the  optimal  policy  obtained  from  above  equations. 
Each  table  reports  the  numerical  results  for  different  cost  and  demand  parameters. 
The  average  costs  presented  in  the  tables  are  taken  over  50  periods  and  over  5,000 
replication.  Figure  3.8  illustrates  the  graphical  results  of  Tables  3.1  ~ 3.4,  respec- 
tively. As  can  be  seen  from  these  tables  and  figures,  the  cost  ratio  of  each  heuristic 
policy  to  the  optimal  policy  varies  with  the  warehouse  capacity.  Except  for  cases 
with  very  tight  warehouse  capacities,  each  heuristic  policy  retains  an  average  cost 
close  to  the  optimal  policy.  We  also  can  observe  that  Heuristic  A has  a relatively 
higher  average  cost  than  Heuristics  B and  C in  most  cases.  Furthermore,  we  can 
see  that  when  the  capacity  increases  and  thus  the  capacity  constraint  becomes  non- 
binding, these  heuristics  yield  near-optimal  expected  costs.  The  latter  result  would 
be  reasonable  because,  as  noted  in  the  previous  section,  each  heuristic  policy  leads 
to  the  unconstrained  optimal  replenishment  policy  when  the  capacity  constraint  is 
not  tight. 

We  now  compare  the  performance  of  each  proposed  heuristic  for  more  general 
cases.  Numerical  trials  are  conducted  for  the  capacity  of  V0  = 100  with  relatively 
low  demand  parameters  compared  to  the  capacity.  In  the  tables  and  figures  below 
we  illustrate  the  results  of  numerical  tests  for  replenishment  schedule  I.  Table  3.5 
presents  the  averages  of  expected  cost  over  all  possible  cost  parameters  for  each 
setting  of  demand  parameters,  while  Figure  3.9  illustrates  the  graphical  results  of 
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Table  3.1:  Heuristic  performance  comparisons  against  the  optimal  policy  for  different 
capacities.  (Ai,A2)  = (30,50),  (hl,h2,pl,p2)  — (2,  2,  2,  2). 


Capacity 

Average  Cost 

A 

Ratio* 

B 

Ratio* 

C 

Ratio* 

Optimal 

20 

288.1 

1.17 

259.9 

1.06 

256.0 

1.04 

246 

40 

254.1 

1.12 

231.5 

1.02 

230.9 

1.02 

227 

60 

230.5 

1.04 

225.9 

1.02 

225.9 

1.02 

223 

80 

221.8 

1.00 

222.7 

1.00 

222.7 

1.00 

222 

100 

221.5 

1.00 

221.5 

1.00 

221.5 

1.00 

221 

120 

221.5 

1.00 

221.5 

1.00 

221.5 

1.00 

221 

Overall 

239.6 

1.05 

230.5 

1.02 

229.8 

1.01 

227 

* Cost/Ratio  to  Optimal 


Table  3.2:  Heuristic  performance  comparisons  against  the  optimal  policy  for  different 
capacities.  (Ai,A2)  = (50,50),  (hl,h2,pl,p2)  — (2,  2,  2,  2). 


Capacity 

Average  Cost 

A 

Ratio* 

B 

Ratio* 

C 

Ratio* 

Optimal 

40 

330.2 

1.19 

288.5 

1.04 

288.5 

1.04 

278.4 

60 

296.7 

1.14 

268.3 

1.03 

269.3 

1.04 

259.5 

80 

273.4 

1.06 

261.2 

1.02 

261.1 

1.02 

257.2 

100 

259.7 

1.01 

257.9 

1.01 

257.9 

1.01 

256.3 

120 

256.3 

1.00 

256.3 

1.00 

256.3 

1.00 

256.3 

Overall 

283.3 

1.08 

266.4 

1.02 

266.6 

1.02 

261.6 

* Cost/Ratio  to  Optimal 


Table  3.3:  Heuristic  performance  comparisons  against  the  optimal  policy  for  different 
capacities.  (Ai,A2)  = (30,10),  (hl,h2,pl,p2)  = (2,  2,  5,  5). 


Capacity 

Average  Cost 

A 

Ratio* 

B 

Ratio* 

C 

Ratio* 

Optimal 

40 

225.3 

1.23 

197.6 

1.08 

197.6 

1.08 

183 

60 

181.3 

1.08 

168.8 

1.01 

169.4 

1.01 

168 

80 

176.5 

1.06 

166.8 

1.00 

166.8 

1.00 

166 

100 

166.3 

1.00 

166.3 

1.00 

166.3 

1.00 

166 

120 

166.3 

1.00 

166.3 

1.00 

166.3 

1.00 

166 

Overall 

183.1 

1.07 

173.2 

1.02 

173.3 

1.02 

170 

* Cost/Ratio  to  Optimal 
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Capacity 


♦ Heuristic.  A 

--»•••  Heuristic.B 
- - Heuristic.C 

Optimal 


(a) 


(b) 


■■  Heuristic_A 
Heuristic_B 
— A — Heuristic_C 
_ - K . - Optimal 


— ♦—Heuristic.  A 
Heurislic.B 
— A — Heuristic.C 
- - * • - Optimal 


Capacity 


(C)  (d) 


Figure  3.8:  Heuristic  performance  comparisons  against  the  optimal  policy  for 

different  capacity,  (a)  (Ai,A2)  = (30,50),  (hl,h2,pl,p2)  = (2,  2,  2, 
2).  (b)  (A1;A2)  = (50,50),  (hl,h2,pl,p2)  = (2,  2,  2,  2).  (c)  (A1;A2) 
= (10,50),  (hl,h2,pl,p2)  = (2,  2,  5,  5).  (d)  (A1;A2)  = (50,10), 

(hi,  h2,pl,p2)  = (2,  2,  20,  20). 
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Table  3.4:  Heuristic  performance  comparisons  against  the  optimal  policy  for  different 
capacities.  (Ai,A2)  = (10,50),  (hl,h2,pl,p2)  = (2,  2,  20,  20). 


Capacity 

Average  Cost 

A 

Ratio* 

B 

Ratio* 

C 

Ratio* 

Optimal 

40 

1815.3 

1.16 

1664.5 

1.07 

1934.5 

1.24 

1558.5 

60 

1594.5 

1.15 

1487.4 

1.07 

1627.9 

1.17 

1389.6 

80 

1460.4 

1.09 

1381.8 

1.04 

1445.2 

1.08 

1333.8 

100 

1372.2 

1.04 

1331.3 

1.01 

1345.8 

1.02 

1321.8 

120 

1323.3 

1.02 

1311.4 

1.01 

1325.6 

1.02 

1301.8 

140 

1299.8 

1.00 

1297.1 

1.00 

1297.1 

1.00 

1296.5 

160 

1295.6 

1.00 

1295.6 

1.00 

1295.6 

1.00 

1295.6 

180 

1295.6 

1.00 

1295.6 

1.00 

1295.6 

1.00 

1295.6 

Overall 

1432.1 

1.06 

1383.1 

1.02 

1445.9 

1.07 

1349.2 

* Cost/Ratio  to  Optimal 


Table  3.5.  Table  3.6  presents  the  average  values  of  expected  cost  over  all  possible  de- 
mand parameters  for  given  specific  cost  parameters,  and  those  results  are  graphically 
represented  in  Figure  3.10.  As  can  be  seen  from  Table  3.5  and  Figure  3.9,  for  different 
demand  rates,  all  three  heuristics  show  almost  identical  overall  performances.  We 
also  can  notice,  from  Table  3.6  and  Figure  3.10,  that  this  result  is  retained  for  differ- 
ent cost  parameters.  This  is  true  because,  in  most  cases,  the  capacity  constraint  is 
not  binding  due  to  low  demand  rates  and  thus  the  problem  becomes  unconstrained. 
Without  the  capacity  constraint,  all  three  heuristics  have  the  identical  replenishment 
policy  which  replenishes  products  up  to  the  optimal  inventory  level  in  the  uncon- 
strained system,  as  mentioned  in  the  previous  section.  This  result  is  consistent  with 
replenishment  schedules  II  and  III  (we  omit  the  results  for  these  schedules).  Since 
it  is  hardly  possible  to  identify  which  heuristic  shows  the  best  performance,  we  next 
consider  the  case  with  the  capacity  being  tight  to  see  if  any  heuristic  outperforms 


the  others. 
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Table  3.5:  Heuristic  performance  for  different  demand  parameters. 


Demand 

Average  Cost 

^1 

^2 

A 

Ratio* 

B 

Ratio* 

C 

10 

10 

1398.3 

1.00 

1398.3 

1.00 

1398.3 

10 

20 

2146.4 

1.00 

2146.4 

1.00 

2146.4 

10 

30 

2892.6 

1.00 

2892.6 

1.00 

2892.6 

20 

10 

2052.3 

1.00 

2052.3 

1.00 

2052.3 

20 

20 

2800.5 

1.00 

2800.5 

1.00 

2800.5 

20 

30 

3581.3 

1.00 

3551.5 

0.99 

3586.6 

30 

10 

2705.5 

1.00 

2705.5 

1.00 

2705.5 

30 

20 

3487.3 

1.00 

3457.4 

0.99 

3479.9 

Overall 

2633.0 

1.00 

2625.6 

1.00 

2632.8 

* Cost /Ratio  to  Heuristic  C 


10  10  10  20  20  20  20  30 

Demands 


- - HeurislicA 

.......  Heuristic_B 

— * — Heuristic_C 


Figure  3.9:  Heuristic  performance  for  different  demand  parameters. 
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Table  3.6:  Heuristic  performance  for  different  cost  parameters. 


Unit  Costs 

Average  Cost 

hi 

h2 

Pi 

P2 

A 

Ratio* 

B 

Ratio* 

C 

2 

2 

2 

2 

2092.6 

1.00 

2092.6 

1.00 

2092.4 

2 

2 

5 

5 

3799.5 

1.00 

3681.4 

0.97 

3787.5 

2 

5 

2 

5 

3573.3 

1.00 

3573.2 

1.00 

3573.2 

5 

2 

5 

2 

3750.9 

1.00 

3750.7 

1.00 

3750.1 

5 

5 

2 

2 

2693.5 

1.00 

2693.5 

1.00 

2693.5 

* Cost /Ratio  to  Heuristic  C 


Figure  3.10:  Heuristic  performance  for  different  cost  parameters. 
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In  new  simulation  runs,  we  vary  the  capacity  while  cost  and  demand  parameters 
are  fixed.  The  choice  of  various  capacities  ranges  from  20  to  400.  This  capacity 
range  seems  to  be  reasonable  for  our  purpose  because  when  the  capacity  is  20,  the 
constraint  is  likely  tight  in  most  cases  with  the  aforementioned  demand  rates,  and 
when  it  is  400,  the  constraint  would  most  likely  be  non-binding.  Henceforth,  we  can 
identify  the  performance  of  each  heuristic  whether  the  capacity  constraint  is  tight  or 
non-binding. 

We  generate  a number  of  problems  from  the  set  of  cost  and  demand  parameters. 
Tables  3.7  ~ 3.11  present  average  costs  for  different  capacities  for  replenishment 
schedule  I.  Tables  3.12  ~ 3.16  report  the  results  for  Schedule  II.  The  results  for 
Schedule  III  are  shown  in  Tables  3.17  ~ 3.21.  In  this  study  we  present  computational 
results  for  five  different  problems  for  each  replenishment  schedule.  The  results  for 
other  problems  retain  a similar  pattern  to  those  in  the  above  tables,  and  thus  we 
omit  them  here. 

In  these  tables,  we  compare  the  performance  of  each  heuristic  to  the  uncon- 
strained optimal  replenishment  policy  with  respect  to  the  expected  cost  over  the 
planning  horizon.  A/U,  B/U  and  C/U  represent  the  cost-ratios  of  heuristics  A,  B 
and  C to  the  unconstrained  optimal  policy,  respectively.  We  also  compare  the  per- 
formance of  heuristics  B and  C to  Heuristic  A,  which  we  consider  as  a base  since  the 
heuristic  is  the  most  conservative  one. 

Note  that  when  there  is  a sufficient  amount  of  slack  in  the  resource  constraint, 
all  heuristics  seem  to  have  identical  performance,  and  thus  the  cost- ratio  of  each 
heuristic  to  the  unconstrained  optimal  policy  approaches  to  one.  This  is  because, 
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when  the  system  becomes  unconstrained,  all  of  them  employ  the  same  replenishment 
policy  that  replenished  products  up  to  the  unconstrained  optimal  replenishment  level. 

As  can  be  seen  from  Tables  3.7  through  3.9  for  Schedule  I,  the  impact  of  tightening 
the  capacity  constraint  is  significant.  When  resources  are  scarce,  Heuristics  B and 
C have  better  performance  than  Heuristic  A (Heuristic  A has  4-11  % higher  average 
cost  than  the  other  two).  This  behavior  is  not  surprising,  especially  since  Heuristic 
A holds  the  smallest  individual  capacity  for  each  item  and  it  is  more  likely  to  retain 
unused  capacity  for  some  items,  while  some  other  items  need  more  storage  to  satisfy 
their  customer  demands,  causing  unnecessary  penalty  costs.  Heuristic  C,  however, 
does  not  carry  this  capacity-inefficiency  incurred  in  Heuristic  A because  it  uses  the 
total  available  capacity  whenever  the  warehouse  replenishes  products.  Heuristic 
B would  confront  the  same  pitfall  as  Heuristic  A because  it  also  adopts  separate 
capacities  for  replenishment,  but  this  flaw  is  lessened  by  the  capacity-adjustment 
characteristic  of  the  heuristic. 

On  the  other  hand,  Tables  3.10  and  3.11  contain  results  for  a special  case  with 
high  penalty  costs  (pi  = 20).  Obviously,  Heuristic  C leads  to  poorer  performance 
when  the  capacity  is  getting  tight.  The  cost  ratio  of  Heuristic  C to  Heuristic  A 
ranges  from  0.98  to  1.18,  depending  on  the  capacity,  while  the  ratio  of  Heuristic  B 
to  Heuristic  A is  as  low  as  0.86.  This  is  most  likely  due  to  the  fact  that,  when  the 
penalty  cost  of  each  product  (pi  = 20)  is  relatively  high  with  respect  to  the  holding 
cost  (hi  = 2),  this  heuristic  is  inclined  to  hold  excessive  inventories  for  products 
being  replenished  in  a period  in  order  to  avoid  the  expensive  penalty  cost,  and  it 
leaves  little  slack  in  the  resource  constraint  in  the  following  replenishment  period, 


81 


causing  high  penalty  costs  for  the  products  to  be  replenished  in  the  next  period. 
Heuristics  A and  B,  however,  are  limited  in  using  any  available  capacity  because  of 
using  separate  capacity  for  each  product,  so  that  both  of  the  heuristics  still  keep 
some  space  for  the  next  replenishment  and  do  not  carry  the  same  degree  of  severity 
as  Heuristic  C does. 


3.6  Conclusions 

In  this  chapter,  we  generalize  the  multi-item,  periodic- review  inventory  model 
described  in  the  previous  chapter,  by  relaxing  the  assumption  of  identical  replenish- 
ment schedules  for  items.  We  first  completely  characterized  the  optimal  replenish- 
ment policy  for  the  two-item,  two-period  case.  For  more  general  cases,  because  of  the 
inability  to  show  convexity  of  the  objective  function,  we  construct  three  heuristics  to 
determine  replenishment  quantities  for  the  case  of  unequal  replenishment  intervals. 

Numerical  testing  of  these  heuristics  suggests  that  they  yield  near  optimal  so- 
lutions for  a small  set  of  problems.  The  performance  of  these  heuristics  has  been 
evaluated  for  a number  of  problems  with  various  cost  and  demand  parameters.  When 
there  is  enough  slack  in  the  resource  constraint,  all  heuristics  demonstrate  identical 
results. 

The  first  heuristic,  denoted  Heuristic  A,  which  uses  a separate  capacity  for  each 
product  for  replenishment,  is  simple  enough  for  practical  use,  but  is  not  able  to  detect 
a scarce  resource  well,  compared  to  the  other  two  heuristics,  in  most  problems.  This 
behavior  is  plausible  because  the  heuristic  is  more  likely  to  retain  unused  capacities 
for  some  products  while  some  others  suffer  from  the  lack  of  the  resource. 
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Heuristic  C adopts  a replenishment  policy  which  replenishes  products  up  to 
the  total  available  capacity  of  the  system  at  the  time  instant  of  replenishment. 
This  heuristic  outperforms  Heuristic  A in  most  cases  because  it  does  not  carry  the 
capacity-inefficiency  confronted  by  the  Heuristic  A , but  leads  to  poorer  performance 
for  the  problems  with  high  penalty  costs.  This  is  because,  in  this  situation,  Heuristic 
C is  too  aggressive  to  replenish  products  in  one  period  so  that  the  warehouse  suffers 
from  the  lack  of  resource  in  the  following  replenishment  period,  causing  high  penalty 
costs  to  the  warehouse. 

On  the  other  hand,  the  second  heuristic,  denoted  Heuristic  B,  adjusts  separate 
capacities  in  each  replenishment  period,  based  on  the  availability  of  any  unused 
resource.  This  characteristic  discriminates  Heuristic  B from  Heuristic  A and  dimin- 
ishes the  capacity-inefficiency  of  Heuristic  A.  Furthermore,  the  capacity-separation 
scheme  of  the  heuristic  leads  to  better  performance  than  Heuristic  C. 
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Table  3.7:  Heuristic  performance  for  different  capacities  with  Schedule  I.  (Ai,A2)  = 
(30,30),  (hl,h2,pl,p2)  = (2,  2,  2,  2). 


Cost/Ratio 

Capacity 

A/U 

B/U 

c/u 

B/A 

C/A 

20 

1.80 

1.60 

1.67 

0.89 

0.93 

40 

1.46 

1.36 

1.41 

0.93 

0.97 

60 

1.21 

1.13 

1.16 

0.94 

0.96 

80 

1.05 

1.05 

1.05 

1.00 

1.00 

100 

1.00 

1.00 

1.00 

1.00 

1.00 

120 

1.00 

1.00 

1.00 

1.00 

1.00 

140 

1.00 

1.00 

1.00 

1.00 

1.00 

Overall 

1.30 

1.23 

1.26 

0.95 

0.97 

Table  3.8:  Heuristic  performance  for  different  capacities  with  Schedule  I.  (Ai,A2)  = 
(30,30),  (hi,  h2,pl,p2)  = (5, 5, 2, 2). 


Cost/Ratio 

Capacity 

A/U 

B/U 

C/U 

B/A 

C/A 

20 

1.39 

1.25 

1.26 

0.90 

0.90 

40 

1.15 

1.08 

1.08 

0.94 

0.94 

60 

1.01 

1.00 

1.00 

0.99 

0.99 

80 

1.00 

1.00 

1.00 

1.00 

1.00 

100 

1.00 

1.00 

1.00 

1.00 

1.00 

120 

1.00 

1.00 

1.00 

1.00 

1.00 

Overall 

1.14 

1.08 

1.08 

0.96 

0.96 

Table  3.9:  Heuristic  performance  for  different  capacities  with  Schedule  I.  (Ai,A2)  = 
(10,50),  (hi,  h2,pl,p2)  = (2,  2,  2,  2). 


Cost/Ratio 

Capacity 

A/U 

B/U 

C/U 

B/A 

C/A 

20 

1.72 

1.64 

1.65 

0.95 

0.96 

40 

1.42 

1.35 

1.36 

0.95 

0.96 

60 

1.19 

1.16 

1.17 

0.97 

0.98 

80 

1.05 

1.03 

1.05 

0.98 

1.00 

100 

1.00 

1.00 

1.00 

1.00 

1.00 

120 

1.00 

1.00 

1.00 

1.00 

1.00 

140 

1.00 

1.00 

1.00 

1.00 

1.00 

Overall 

1.28 

1.24 

1.24 

0.97 

0.98 
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Table  3.10:  Heuristic  performance  for  different  capacities  with  Schedule  I.  (Ai,  A2)  = 
(10,50),  (hi,  h2,pl,p2)  = (2,  2,  20,  20). 


Cost/Ratio 

Capacity 

A/U 

B/U 

c/u 

B/A 

C/A 

20 

5.18 

4.90 

5.69 

0.94 

1.10 

40 

4.21 

3.94 

4.71 

0.94 

1.12 

60 

3.37 

3.16 

3.75 

0.94 

1.11 

80 

2.67 

2.51 

2.95 

0.94 

1.10 

100 

2.12 

2.02 

2.30 

0.95 

1.08 

120 

1.72 

1.65 

1.82 

0.96 

1.06 

140 

1.42 

1.38 

1.49 

0.97 

1.04 

160 

1.22 

1.21 

1.26 

0.99 

1.03 

180 

1.09 

1.09 

1.12 

0.99 

1.02 

200 

1.02 

1.02 

1.04 

1.00 

1.02 

220 

1.00 

1.00 

1.00 

1.00 

1.00 

240 

1.00 

1.00 

1.00 

1.00 

1.00 

260 

1.00 

1.00 

1.00 

1.00 

1.00 

Overall 

2.17 

2.07 

2.34 

0.97 

1.06 

Table  3.11:  Heuristic  performance  for  different  capacities  with  Schedule  I.  (Ai,  A2)  = 
(30,30),  (hl,h2,pl,p2)  = (2,  2,  20,  20). 


Cost/Ratio 

Capacity 

A/U 

B/U 

C/U 

B/A 

C/A 

Capacity 

A 

Ratio* 

B 

Ratio* 

C 

20 

5.48 

4.85 

5.39 

0.89 

0.98 

40 

4.42 

3.79 

4.55 

0.86 

1.03 

60 

3.50 

3.01 

3.81 

0.86 

1.09 

80 

2.75 

2.40 

3.14 

0.87 

1.14 

O 

O 

r— 1 

2.17 

1.93 

2.56 

0.89 

1.18 

120 

1.74 

1.58 

2.04 

0.91 

1.17 

140 

1.44 

1.33 

1.58 

0.92 

1.10 

160 

1.23 

1.16 

1.30 

0.95 

1.06 

180 

1.09 

1.07 

1.13 

0.98 

1.04 

200 

1.02 

1.02 

1.04 

1.00 

1.02 

220 

1.00 

1.00 

1.00 

1.00 

1.00 

240 

1.00 

1.00 

1.00 

1.00 

1.00 

260 

1.00 

1.00 

1.00 

1.00 

1.00 

Overall 

2.24 

2.01 

2.38 

0.93 

1.07 
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Table  3.12:  Heuristic  performance  for  different  capacities  with  Schedule  II.  (Ai,A2) 
= (30,30),  (hl,h2,pl,p2)  = (2,  2,  2,  2). 


Cost/Ratio 

Capacity 

A/U 

B/U 

c/u 

B/A 

C/A 

40 

1.97 

1.78 

1.72 

0.90 

0.87 

60 

1.65 

1.57 

1.57 

0.95 

0.95 

80 

1.56 

1.46 

1.46 

0.94 

0.94 

100 

1.52 

1.46 

1.46 

0.96 

0.96 

120 

1.11 

1.07 

1.07 

0.97 

0.97 

140 

1.07 

1.04 

1.04 

0.97 

0.97 

160 

1.00 

1.00 

1.00 

1.00 

1.00 

180 

1.00 

1.00 

1.00 

1.00 

1.00 

Overall 

1.41 

1.34 

1.33 

0.96 

0.95 

Table  3.13:  Heuristic  performance  for  different  capacities  with  Schedule  II.  (Ai,A2) 
= (10,50),  (hl,h2,pl,p2)  = (5,  5,  2,  2). 


Cost/Ratio 

Capacity 

A/U 

B/U 

C/U 

B/A 

C/A 

40 

1.77 

1.71 

1.68 

0.96 

0.95 

60 

1.61 

1.53 

1.51 

0.96 

0.94 

80 

1.43 

1.36 

1.36 

0.95 

0.95 

100 

1.34 

1.31 

1.31 

0.98 

0.98 

120 

1.33 

1.31 

1.31 

0.99 

0.99 

140 

1.02 

1.01 

1.01 

0.99 

0.99 

160 

1.00 

1.00 

1.00 

1.00 

1.00 

Overall 

1.36 

1.32 

1.31 

0.98 

0.97 

Table  3.14:  Heuristic  performance  for  different  capacities  with  Schedule  II.  (Ai,A2) 
= (50,30),  (hl,h2,pl,p2)  = (2,  2,  2,  2). 


Cost/Ratio 

Capacity 

A/U 

B/U 

C/U 

B/A 

C/A 

60 

1.77 

1.59 

1.57 

0.90 

0.88 

80 

1.54 

1.44 

1.45 

0.94 

0.94 

100 

1.45 

1.36 

1.36 

0.94 

0.94 

120 

1.24 

1.19 

1.19 

0.96 

0.97 

140 

1.11 

1.09 

1.09 

0.98 

0.98 

160 

1.04 

1.04 

1.04 

1.00 

1.00 

180 

1.00 

1.00 

1.00 

1.00 

1.00 

Overall 

1.31 

1.24 

1.24 

0.96 

0.96 

86 


Table  3.15:  Heuristic  performance  for  different  capacities  with  Schedule  II.  (Ai,A2) 
= (30,10),  (hi,  h2,pl,p2)  = (2,  2,  20,  20). 


Cost/Ratio 

Capacity 

A/U 

B/U 

c/u 

B/A 

C/A 

60 

3.47 

3.13 

4.65 

0.90 

1.34 

80 

2.97 

2.73 

3.46 

0.92 

1.16 

100 

2.73 

2.57 

2.89 

0.94 

1.06 

120 

1.33 

1.30 

1.49 

0.97 

1.12 

140 

1.11 

1.13 

1.24 

1.02 

1.12 

160 

1.02 

1.06 

1.07 

1.04 

1.05 

180 

1.01 

1.01 

1.01 

1.01 

1.01 

200 

1.00 

1.00 

1.00 

1.00 

1.00 

220 

1.00 

1.00 

1.00 

1.00 

1.00 

Overall 

1.83 

1.74 

2.10 

0.98 

1.11 

Table  3.16:  Heuristic  performance  for  different  capacities  with  Schedule  II.  (Ai,A2) 
= (50,30),  (hi,  h2,pl,p2)  = (2,  2,  20,  20). 


Cost/Ratio 

Capacity 

A/U 

B/U 

C/U 

B/A 

C/A 

60 

3.27 

2.91 

3.65 

0.89 

1.11 

80 

2.61 

2.33 

2.88 

0.89 

1.10 

100 

2.11 

1.91 

2.24 

0.91 

1.06 

120 

1.72 

1.58 

1.79 

0.92 

1.05 

140 

1.47 

1.40 

1.49 

0.95 

1.02 

160 

1.41 

1.37 

1.36 

0.97 

0.96 

180 

1.10 

1.10 

1.17 

1.00 

1.06 

200 

1.03 

1.04 

1.08 

1.02 

1.05 

220 

1.01 

1.03 

1.02 

1.02 

1.01 

240 

1.00 

1.01 

1.01 

1.01 

1.01 

260 

1.00 

1.00 

1.00 

1.00 

1.00 

280 

1.00 

1.00 

1.00 

1.00 

1.00 

Overall 

1.61 

1.52 

1.70 

0.96 

1.04 
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Table  3.17: 


Heuristic  performance  for  different  capacities  with  Schedule  III. 
(Ai,A2,A3)  = (30,10,30),  (hl,h2,h3,pl,p2,p3)  = (2,  2,  2,  2,  2,  2). 


Cost/Ratio 

Capacity 

A/U 

B/U 

c/u 

B/A 

C/A 

50 

1.55 

1.34 

1.41 

0.87 

0.91 

70 

1.41 

1.29 

1.28 

0.91 

0.91 

90 

1.19 

1.07 

1.14 

0.90 

0.96 

110 

1.08 

1.02 

1.05 

0.94 

0.97 

130 

1.02 

1.00 

1.01 

0.98 

0.99 

150 

1.00 

1.00 

1.00 

1.00 

1.00 

170 

1.00 

1.00 

1.00 

1.00 

1.00 

190 

1.00 

1.00 

1.00 

1.00 

1.00 

Overall 

1.18 

1.10 

1.13 

0.94 

0.96 

Table  3.18: 


Heuristic  performance  for  different  capacities  with  Schedule  III. 
(Ai,A2,A3)  = (10,10,50),  (hl,h2,h3,pl,p2,p3)  = (2,  2,  2,  2,  2,  2). 


Cost/Ratio 

Capacity 

A/U 

B/U 

C/U 

B/A 

C/A 

50 

1.61 

1.42 

1.49 

0.88 

0.93 

70 

1.41 

1.26 

1.31 

0.89 

0.93 

90 

1.26 

1.14 

1.17 

0.91 

0.93 

110 

1.15 

1.12 

1.12 

0.97 

0.97 

130 

1.13 

1.11 

1.11 

0.99 

0.99 

150 

1.04 

1.02 

1.02 

0.99 

0.99 

170 

1.00 

1.00 

1.00 

1.00 

1.00 

190 

1.00 

1.00 

1.00 

1.00 

1.00 

Overall 

1.23 

1.15 

1.17 

0.95 

0.96 

Table  3.19: 


Heuristic  performance  for  different  capacities  with  Schedule  III. 
(Ai,A2,A3)  = (10,10,30),  (hl,h2,h3,pl,p2,p3)  = (2,  2,  2,  5,  5,  5). 


Cost/Ratio 

Capacity 

A/U 

B/U 

C/U 

B/A 

C/A 

30 

2.43 

2.03 

2.19 

0.84 

0.90 

50 

2.01 

1.74 

1.80 

0.87 

0.90 

70 

1.80 

1.64 

1.63 

0.91 

0.91 

90 

1.69 

1.59 

1.55 

0.94 

0.92 

110 

1.21 

1.09 

1.23 

0.90 

1.02 

130 

1.08 

1.02 

1.12 

0.95 

1.04 

150 

1.02 

1.01 

1.02 

0.99 

1.00 

170 

1.00 

1.00 

1.00 

1.00 

1.00 

190 

1.00 

1.00 

1.00 

1.00 

1.00 

Overall 

1.53 

1.39 

1.44 

0.92 

0.96 
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Table  3.20: 


Heuristic  performance  for  different  capacities  with  Schedule  III. 
(Ai,A2,A3)  = (10,10,10),  (hl,h2,h3,pl,p2,p3)  = (2,  2,  2,  20,  20,  20). 


Cost/Ratio 

Capacity 

A/U 

B/U 

C/U 

B/A 

C/A 

30 

4.26 

3.40 

5.38 

0.80 

1.26 

50 

3.02 

2.41 

4.01 

0.80 

1.33 

70 

2.14 

1.74 

2.99 

0.81 

1.40 

90 

1.58 

1.33 

2.67 

0.84 

1.69 

110 

1.21 

1.10 

1.59 

0.91 

1.31 

130 

1.03 

1.01 

1.03 

0.98 

1.00 

150 

1.00 

1.00 

1.00 

1.00 

1.00 

Overall 

2.04 

1.71 

2.67 

0.88 

1.28 

Table  3.21: 


Heuristic  performance  for  different  capacities  with  Schedule  III. 
(Ai,A2,A3)  = (30,30,30),  {hl,h2,h3,pl,p2,p3)  = (2,  2,  2,  20,  20,  20). 


Cost/Ratio 

Capacity 

A/U 

B/U 

C/U 

B/A 

C/A 

50 

5.28 

4.29 

5.23 

0.81 

0.99 

70 

4.75 

3.82 

4.71 

0.80 

0.99 

90 

4.26 

3.40 

4.22 

0.80 

0.99 

110 

3.83 

3.06 

3.72 

0.80 

0.97 

130 

3.37 

2.79 

3.35 

0.83 

0.99 

150 

3.04 

2.58 

3.00 

0.85 

0.99 

170 

2.70 

2.16 

3.43 

0.80 

1.27 

190 

2.39 

1.93 

3.21 

0.81 

1.34 

210 

2.12 

1.73 

3.03 

0.82 

1.43 

230 

1.91 

1.58 

2.89 

0.82 

1.51 

250 

1.71 

1.44 

2.78 

0.84 

1.62 

270 

1.56 

1.33 

2.68 

0.85 

1.72 

290 

1.41 

1.23 

2.39 

0.87 

1.69 

310 

1.30 

1.16 

1.96 

0.89 

1.50 

330 

1.22 

1.10 

1.58 

0.91 

1.30 

350 

1.14 

1.06 

1.30 

0.93 

1.14 

370 

1.09 

1.03 

1.11 

0.95 

1.02 

390 

1.04 

1.01 

1.03 

0.97 

0.99 

410 

1.01 

1.00 

1.00 

0.99 

0.99 

430 

1.00 

1.00 

1.00 

1.00 

1.00 

450 

1.00 

1.00 

1.00 

1.00 

1.00 

Overall 

2.31 

1.94 

2.68 

0.87 

1.22 

CHAPTER  4 

CUSTOMER  ORDER  LEAD-TIME  MODELS 
4.1  Introduction 

Uncertainties  in  manufacturing  environments  play  important  roles  in  determin- 
ing planned  customer  order  leadtimes.  These  uncertainties  are  mostly  due  to  com- 
plicated production  processes,  random  yields,  and  high  quality  requirements.  In- 
corporation of  these  factors  into  the  selection  of  planned  leadtimes  represents  an 
important  step  toward  increasing  the  robustness  of  manufacturing  planning.  The 
leadtime  uncertainty  itself  affects  many  aspects  of  costs  and  control.  It  is  particu- 
larly problematic  in  queueing  systems  because  tardiness  in  earlier  orders  may  delay 
subsequent  orders. 

Most  of  the  due-date-setting  models  assume  that  the  due  dates  for  individual  or- 
ders are  set  entirely  exogenously.  However,  in  certain  practical  contexts,  each  order 
arrives  with  a due  date,  indicating  some  future  time  when  the  customer  wishes  to 
receive  the  goods  ordered  and,  in  most  cases,  due  date  setting  is  negotiable  and  is 
the  responsibility  of  the  marketing  personnel  of  the  firm.  For  the  marketing  group, 
who  have  knowledge  of  customer’s  wishes,  it  is  important  to  know  their  manufactur- 
ing personnel’s  perspective  on  the  order  before  negotiating  the  due  date  with  their 
customer.  That  is,  they  would  like  to  know  the  optimal  customer  order  leadtime 
to  quote  for  the  order.  By  customer  order  leadtimes,  we  mean  the  time  from  a 
customer’s  order  until  the  due  date. 
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This  study  is  stimulated  by  Yano’s  work  on  PLTs  in  serial  production  systems 
(1987a).  Yano  developed  a solution  approach  to  solve  the  two-stage  problem  and 
applied  it  to  iV-stage  serial  systems.  The  limitation  of  their  approach  was  that 
they  did  not  address  the  effects  of  queueing  or  capacity  explicitly  in  the  model.  For 
instance,  Yano  assumes  that  an  order  that  has  completed  processing  in  one  stage  may 
be  held  until  its  due  date  (i.e.  no  early  delivery  to  customers  is  allowed).  However, 
it  is  also  assumed  that  once  that  order  arrives  at  the  next  stage,  it  is  processed 
immediately  without  being  delayed  by  preceded  orders.  This  assumption  greatly 
simplifies  their  solution  approach  for  determining  PLTs  in  serial  production  systems. 
The  assumption  used  by  Yano  is  not  pragmatic  in  many  situations.  Literature  such 
as  Duenyas  (1995)  and  Duenyas  and  Hopp  (1995)  considered  queueing  effects  in 
determining  PLTs.  In  this  study,  we  consider  the  aforementioned  delay  as  well  as  the 
possible  holding  cost  incurred  due  to  the  prohibition  of  early  delivery  to  customers 
in  single-stage  systems. 

This  study  addresses  single-stage  make-to-order  (MTO)  production  systems.  Our 
main  objective  is  to  determine  the  optimal  planned  customer  leadtimes  while  mini- 
mizing the  sum  of  expected  inventory-holding  costs  and  expected  penalty  costs  for 
exceeding  the  quoted  due  date,  so  as  to  quote  planned  leadtimes  to  customers  at 
the  time  of  order  arrival.  Since  the  area  of  modeling  leadtime  setting  where  demand 
is  sensitive  to  quoted  leadtimes  is  relatively  unexposed,  it  would  not  be  a critical 
shortcoming  to  assume  that  the  leadtimes  quoted  to  customers  are  independent  of 
each  other.  Within  this  scope,  we  focus  on  modeling  the  appropriate  expected  costs 
under  a variety  modeling  assumptions  and  characterizing  the  optimal  policies. 
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Our  approach  to  solving  the  problem  is  to  derive  the  distribution  of  actual  com- 
pletion time  of  the  process  for  an  individual  order  and  to  compare  this  to  the  cor- 
responding quoted  due  dates  to  obtain  the  expected  total  costs.  We  then  minimize 
costs  with  respect  to  the  decision  variable(s),  which,  in  this  case,  are  PLT(s).  We 
show  that,  in  special  cases,  the  single  stage  model  is  equivalent  to  a simple  newsboy 
problem  as  Weeks  (1981)  pointed  out. 


4.2  Model  Description  and  Formulation 
Consider  the  problem  of  determining  optimal  planned  leadtimes  in  single  stage 
make-to-order  production  systems  where  individual  orders  wait  in  a G/G/l  queue 
until  preceded  orders  complete  their  processes.  The  production  system  consists  of 
one  reliable  machine  and  produces  a single  item  (see  Figure  4.1).  Once  the  process 
for  an  order  has  been  completed,  the  finished  goods  are  held  in  inventory  until  the 
due  date.  If  a job  has  already  been  tardy,  the  products  are  dispatched  immediately 
(see  Figure  4.2).  Therefore,  there  are  two  possible  delays,  one  by  holding  products 
until  the  due  date  and  the  other  by  waiting  until  the  previous  order  is  completed. 
The  order  size  is  a random  variable. 

Let  k be  the  index  of  newly  arrived  order.  Denote 
xk  : planned  customer  leadtime  of  the  /cth  order, 

R : random  variable  representing  yield  rates  with  a density  <?i(r),  R > 0, 

Qk  : the  size  of  order  k,  where  Qk  ~ iid  a random  variable  with  a density 

92(q),  Qk  = 1,2,---, 

tk  : arrival  time  of  order  k,  a constant, 

r:  processing  time  of  unit  product  with  a density  <73  (r), 

sk  : actual  process  completion  time  of  the  kth  order, 

h : holding  cost  per  unit  product  and  per  unit  time, 

p : penalty  cost  for  late  delivery  per  unit  product  and  per  unit  time. 
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Figure  4.1:  One  machine,  single  part-type  system. 

Then  the  due  date  of  order  k,  Dk,  becomes  Dk  — xk  +tk  and  the  total  processing 
time  of  the  kth  order,  rk,  is  rk  = ^r.  The  term  ^ represents  the  actual  amount 
of  finished  goods  to  produce  to  meet  customer  demands.  The  arrival  time  of  that 
order,  tk,  is  known  and  constant.  The  actual  process  starting  time  of  order  k is 
max  {tk,  S/t-i},  where  s*-i  is  the  job  completion  time  of  previous  ( k — l)th  order  (see 
Figure  4.2).  This  is  true  because  if  there  is  no  order  in  the  system,  the  newly  arrived 
one  will  be  served  immediately.  Otherwise,  it  must  wait  until  the  preceding  order, 
order  ( k — 1),  is  completed.  Then,  the  actual  process  completion  time  of  order  k,  sk, 
can  be  represented  as  follows: 

sk  = m&x{tk,sk_l}  + Tk  (4.1) 

Since  an  individual  customer  will  not  take  a delivery  before  the  due  date,  a product 
holding  cost  is  incurred  for  early  completion.  The  cost  to  customers  of  early  deliveries 
are  widely  appreciated,  partly  due  to  the  JIT  movement.  On  the  other  hand,  a tardy 
job  will  be  released  immediately  after  completion  with  a penalty  cost.  Under  the 
above  assumptions,  we  minimize  the  expected  cost  function  with  respect  to  the 
customer  planned  leadtime  of  order  k,  xk, 

min  ip(xk)  = E[hQk(Dk  - sk)+  +pQk(sk  - Dk)+ ] (4.2) 

Xk 
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where 


Figure  4.2:  Configuration  of  process. 


a+  = max  {a,  0} 

E[hQk(Dk  - sk)+) 

= E[hQk 

= E[hQk 

E[hQk{sk  - Dky\ 

= E\pQk 

{tk  + *k  ~ Sfc)+] 

{tk  + xk  - m&x{tk,  sk-i}  - rk)+] 

(Sfc  ~ tk  — Xk)  ] 


- E\pQk(max{tk,sk-i} + Tk-tk- xk)+] 


The  expectation  is  taken  with  respect  to  the  random  variables  R,Qk,r  and  sk.  The 
first  term  in  the  right  hand  side  in  (4.2)  denotes  the  expected  holding  cost  due  to 
early  production  and  the  second  one  represents  the  expected  cost  for  tardiness.  To 
solve  this  nonlinear  optimization  problem,  we  first  derive  the  distribution  of  sk. 


4.3  Distribution  of  sk 

Since  sk  is  a function  of  max  { tk , we  consider  two  possible  cases  separately. 

The  first  case  is  such  that  the  system  is  empty  when  order  k arrives,  or  tk  > sk_i. 
In  this  case,  the  process  completion  time  of  the  ( k — l)th  order  is  a known  constant 
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and  Sk  is  just  sum  of  a random  variable,  Tk,  and  a constant  tk. 

Theorem  4.1  Sk  has  the  same  distribution  as  Tk  with  a location  parameter  tk,  where 
tk  > rk. 

Proof  : From  (4.1),  we  get 

sk  = max.{tk,sk-i}  +rk  = tk  + rfc  □ 


Since  tk  is  a known  constant  and  the  distribution  of  Tk  is  given,  we  can  easily  derive 
the  distribution  of  sk. 

Now  we  can  simplify  the  expected  cost  function  (4.2)  as  follows: 

ip(xk)  = E [hQk(tk  + xk  - max{tfc,sfc_i}  - rfc)+] 

+E  [pQk{m&x{tk,  sk-i}  + Tk  -tk-  zfc)+] 

= E [hQk{xk  - Tk)+]  + E [pQk(rk  - xk)+]  (4.3) 

The  second  case  is  such  that  tk  is  less  than  or  equal  to  sk~ i-  In  other  words,  there 
is  at  least  one  customer  order  in  the  system  when  order  k arrives.  Without  lose  of 
generality,  let’s  assume  that  l orders,  including  the  newly  arrive  order  k,  are  in  the 
system  at  time  tk.  Then,  sk  is  a function  of  the  total  processing  time  of  the  last  k — 1 
orders,  the  remaining  processing  time  of  the  order  in  process  at  time  tk  and  a known 
constant,  sk-i,  which  is  the  process  completion  time  of  the  latest  order  to  finish  at  tk- 
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Theorem  4.2  The  distribution  of  Sk  is  a convolution  of  a constant  and  random 
variables  whose  distributions  are  known 

Proof : Since  there  are  l orders  in  the  system,  we  get  tk  < sk,  tk  < sk- i,  ... , tk  < 
Sfc-i+i  and  tk  > sk-i.  The  first  l inequalities  result  from  the  fact  that  the  last  l orders 
in  the  system  have  not  been  served  yet  at  time  tk.  Thus,  the  job  completion  time 
of  those  orders  must  greater  than  tk.  The  last  inequality  is  obvious  because  the 
(k  — J)th  order  has  already  been  processed  at  time  tk  (the  (k  — l + l)th  order  is  cur- 
rently on  processing)  and,  thus,  is  a known  constant  at  time  tk.  We  also  know 
that  ffc-j+i  < tfc-z+2  < • • • < tjfe-i  < tk.  Combining  these  two  sets  of  inequalities,  we 
get 

tk-i+ 1 < sjt_z,  4-z+2  < sk-i+ 1)  •••  tk  < sfc-i  (4-4) 

Applying  (4.4)  to  (4.1),  we  have 

sk  = max{tfc,sfe_i}  + rfc  = sfc-i  + rk 
sk~  i = max{tfc_i,sfc_2}  + Tfc_i  = s*;_2  + 'Tfc_i 


Sfc-Z+l 


ax  ft 


>]p  — 1 -L  1 « & hr  . 


T + 


^ 2 , 

sk  = Sk-l  + > QTk-i  + Tjfc-Z+l 


(4.5) 


where  r'k_l+l  — U + rk_i+i  — tk  is  the  remaining  process  time  of  order  (k  — l + l) 


at  time  tk,  U is  the  actual  process  starting  time  of  order  (k  — l + 1).  Note  that 
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the  job  completion  time  of  order  k,  sk,  is  a function  of  the  job  completion  time  of 


the  (k  - i)th  order,  sk-i,  and  the  projected  total  processing  time  of  orders  waiting 
in  queue  at  time  tk.  Now,  we  can  derive  the  distribution  of  sk  from  (4.5)  since  all 
random  variables  included  in  the  right  hand  side  in  (4.5)  have  known  distributions 
and  is  a known  constant.  □ 

Denoting  the  density  of  sk  be  g4(sk)  we  can  rewrite  the  expected  cost  function 
as  follows: 


be  solved  separately  for  each  of  aforementioned  cases.  We  establish  a closed-form 
solution  for  each  case  and  investigate  the  effects  of  changing  the  parameters  in  our 
model  on  the  optimal  planned  customer  leadtime. 

4.4.1  Case  I:  tk  > sk_ i 

In  case  that  there  are  no  orders  waiting  or  on  processing  in  the  system  when  a 
new  order  arrives  at  time  tk,  we  get,  from  (4.3), 


ip(xk)  = E [hQt(tt  +xt  - st)+]  +E  [pQk(st  -tk  -x*)+] 


(4.6) 


4.4  Solution  Approach 


Once  we  derive  distribution  of  sk,  the  nonlinear  optimization  problem  (4.2)  can 


ip(xk)  = E [hQk(xk  - rjfc)+]  + E [pQk(rk  - xk)+ ] 
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= E 
= h 


hQk  [xk-T 


Qk 

R 


+ E 


n ( Qk 

pQk  It—  - xk 


[ 9i{r)  [ qk92{qk ) f (xk  ~ — ) gz{r)dTdqkdr 
Jr  Jqk  Jt= 0 V T J 

+P  f 9i(r)  [ qk92(qk)  f (—  - xk)  9z{r)dTdqkdr 
Jr  Jqk  Jt= 0 V r ' 


= fc  fSd D [ M r'n  (x,  - 

Jr  Jqk  Jt= 0 X r 

+P  f 9i(r)  [ 9k92(qk)  [ C^r  ~ xk)  9z{r)dTdqkdr  (4.7) 

Jr  Jqk  J T=rxk/qk  ' r 

and  we  want  to  minimize  (4.7)  with  respect  to  the  planned  customer  leadtime,  xk- 
Assuming  ij}(xk)  is  twice  differentiable,  we  can  then  prove  the  following: 


g3(T)drdqkdr 


Proposition  4.1  ip(xk)  is  convex  in  Xk  > 0 

Proof  : Differentiating  (4.7)  with  respect  to  xk,  we  obtain  the  following  optimality 
conditions: 


dx 


mxk)  = h l 91  w L £ -t  [IS  (xh~^  9s(t)* 

+P  9i(r)  qk92{qk)w—  / ~ xkJ  g3{r)dT  dqkdr 

Jr  J qk  0Xk  \-J r=rxk/qk  V r 


dqkdr 


r r rTxk/qk 

= h gi{r)  qk92{qk ) / g3{r)dTdqkdr 

Jr  J qk  J T— 0 

+P  [ 9i(r)  [ qk92(qk ) [ -g3{T)drdqkdr 

Jr  J ql  J T=rxk/Qk 


r r \ rTXk/vk 

= 9i(r)  qk92{qk)\h  l g3(r)dT  - p gz{r)dn 

Jr  Jqk  [ ^ T=°  J T—rxk/qk 


dqkdr 
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= 

f 9i(r) 

1 qk92(qk) 

rrxk/qk  roc 

(h  + p)  g3(r)dr  — pi  g3(r)dr 

J 

r «/ 

<ik 

1 

O 

II 

h 

o 

II 

h 

= 1 

f 9i{r) 

j Qkg2(qk) 

rrxk/qk 

(h  + p)  g3(r)dr  — p\  dqkdr 

J 

r J 

Qk 

Jt= 0 

rrxk/Qk 

= (h+p)  J 

9i(r)  / qkg2(qk)  / g3(r)drdqkdr  - pE[Qk] 

r Jqk  J t= 0 

dqkdr 


= (h  + p)  g i(r)  / qkG3(rxk/qk)g2(qk)dqkdr -pE[Qk]  = 0 (4.8) 

7r  Jqk 


d2ip{xk) 
dxk 2 


92{qk)dqkdr 


= {h  + p)  Ir91^  Iq  qk  \^kG^rXk/qk^ 

= (h  + p)  gi{r)  qk(r/qk)g3(rxk/qk)g2(qk)dqkdr 
Jt  Jqk 

= (h  + p)  / rgi(r)  / g3(rxk/qk)g2(qk)dqkdr 
Jr  J au 


(4.9) 


where  G3(-)  is  a cumulative  distribution  function  of  g3(-)-  Note  that  the  right  hand 
side  in  (4.9)  is  nonnegative  because  all  terms  inside  integrations  are  nonnegative. 
Thus,  the  expected  cost  function  is  convex  and  the  solution  satisfying  the  first-order 
condition  provides  an  unique  minimum.  □ 

Toward  this  end,  we  now  solve  (4.8)  under  the  assumption  that  the  order  size  is 
uniformly  distributed  where  0 < qk  < N,  the  random  yield  ranges  uniformly  between 
zero  and  one,  and  the  unit  processing  time  has  a density  of  uniform,  17(0,  T).  Under 
the  above  assumptions,  we  state  the  following: 


Proposition  4.2  If  gi(r)  = 1,  g2(qk)  - 1/N  and  g3(r)  = l/T, 
then  the  optimal  solution  to  (4.7)  is 

x*  = ThTp)E[Qk]i~2T^  where  E ^ = Y 
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Proof  : Substituting  gi(r),g2(qk)  and  </3(t)  into  (4.8)  and  solving  for  xk  gives 


Therefore,  in  case  that  there  is  no  order  in  the  system  when  the  kth.  order  arrives 
at  time  tk,  the  optimal  solution,  x*k,  is  a function  of  the  maximum  processing  time, 
average  order  size  as  well  as  the  holding  cost  and  the  tardiness  cost.  This  solution 
along  with  N = 1 and  T = 1 establishes  its  equivalence  with  a simple  newsboy 
problem  as  Weeks  (1981)  showed.  That  is,  in  case  that  no  orders  are  waiting  or 
being  processed  in  the  system,  our  single  stage  problem  is  equivalent  to  a simple 
newsboy  problem.  Next,  we  consider  the  case  that  at  least  one  order  is  in  the  system 
when  a new  order  arrives  at  time  tk. 

4.4.2  Case  II:  tk  < sk_ i 

This  is  the  case  that  when  a new  customer  order  arrives  at  time  tk,  at  least  one 
preceded  order  is  in  the  system.  Referring  to  (4.6),  we  have 
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= hf  qkg2(qk)  [ (Ak  - sk)  gi(sk)dskdqk 
Jqk  Jsk=0 

+P  qkg2(qk)  (sk  - Ak)  g4(sk)dskdqk 

J Qk  J Sk—Au 


(4.10) 


where  Ak  = tk  + xk  and  g4  is  the  density  function  of  sk. 

We  want  to  minimize  (4.10)  with  respect  to  the  planned  customer  leadtime,  xk. 
Again,  assuming  ip(xk)  is  twice  differentiable,  we  can  prove  the  convexity  of  ip(xk) 
as  follows: 


Proposition  4.3  ip(xk)  is  convex  in  xk  > 0 
Proof : Prom  (4.10),  it  follows  that 


dip{xk) 

dxk 


— h f qkg2{qk)  n [ [ (Ak  sk)  g4(sk)dsk  dqk 

Jqk  Usk=  0 J 

+P  / qkg2{qk )«—  / (sk  - Ak)  g4(sk)dsk  dqk 

Jqk  OXk  lJsk=Ak 

=h  qkg2(qk)  / g4(sk)dsk  dqk 

Jqk  lJsk=0 

+p  / qkg2(qk ) / -g4(sk)dsk  \ dqk 

Jqk  Usk=Ak  J 

r r Mk  r°° 

= / Qkg2(qk)\h  94(sk)dsk-p  / g4(sk)dsk  dqk 

Jqk  L J sk=  0 J sk=Ak 

r ["  rAk  roo 

= / Qk92{qk)\{h  + p)  I g4{sk)dsk-p  g4{sk)dsk  dqk 

Jqk  L J sk=0  Jsk=0 


= [ qk92{qk ) (h  + p)  f g4(sk)dsk  — p dqk 

Jqk  L Jsk=  0 

= / qk[{h  + p)G4(Ak)  - p\g2(qk)dqk 

J Qk 


(4.11) 


where  G^-)  is  a cumulative  distribution  function  of  sk.  Equating  (4.11)  to  zero,  the 
optimal  solution  x*k  satisfies  the  following  first  order  condition 
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(h  + p)  / qkG4(tk  + x*k)g2(qk)dqk  - pE[Qk } = 0 


(4.12) 


and  the  second  derivative  is 


d2ip(xk) 
dxk 2 


Qk 


which  is  nonnegative  since  all  terms  inside  the  integrations  are  nonnegative.  Thus, 
the  expected  cost  function  (4.10)  is  convex  and  the  solution  satisfying  the  first-order 
condition  (4.12)  provides  the  unique  optimal.  □ 

Unfortunately,  the  solution  of  (4.12)  is,  in  many  cases,  not  tractable  analytically 
and  we  thus  employ  an  numerical  approximation  to  retrieve  the  solution. 


In  this  section,  we  study  our  model  numerically  to  investigate  the  effect  of  chang- 
ing values  of  cost  parameters  and  average  processing  time  on  the  optimal  expected 
cost.  For  simplicity,  we  assume  that  each  customer  orders  a single  unit  ( Qk  = 1) 
and  that  there  is  no  yield  variability.  We  also  assume  that  there  are  only  two  orders 
in  the  system  when  the  kth  order  arrives  at  time  tk.  The  actual  process  completion 
time  of  order  k is  then,  from  (4.5), 


4.5  Numerical  Results 


Sk  — ■Sfc_2  + Tfc  + rfc_1 


Sk-2  + Tk  + Tfc_i  +U  — tk 
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where  U is  the  actual  process  starting  time  of  order  (k  — 1).  If  the  actual  processing 
time  of  individual  orders  is  exponentially  distributed  with  a mean  1/A,  the  density 
function  of  sk  is 


c = sfc_2  + U - tk.  Referring  to  (4.12),  we  get  the  following  first-order  condition: 


Since  Sfc_2,  U and  tk  are  all  known  constants,  we  set  s*_2  = 5.0,  U = 3.0  and 
tk  = 4.0.  The  cost  parameters  considered  for  this  problem  axe  as  follows:  unit 
holding  costs  range  h — 1,2,..., 10  and  unit  penalty  costs  are  p = 1, 2, ...,  10.  We  also 
vary  mean  processing  times,  A,  between  one  time  unit  and  fifty. 

Numerical  trials  are  conducted  for  a number  of  different  parameter  settings.  For 
each  parameter  setting,  we  use  simulation  to  estimate  the  optimal  PLT,  x*k  and  ex- 
pected cost,  from  (4.11).  Figures  4.3  through  4.5  contain  examples  of  our 

findings. 

1.  Effect  of  average  processing  time  on  optimal  PLT  and  expected  cost.  Figure  4.3 
presents  the  relationship  between  the  average  processing  time  and  the  optimal  PLT, 
and  between  the  average  processing  time  and  the  optimal  cost,  given  that  the  value 
of  cost  parameters  are  fixed.  As  we  can  see,  both  optimal  PLT  and  expected  cost 
increase  with  an  increased  processing  time.  This  result  is  plausible  because  a longer 
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Figure  4.3:  Optimal  cost  and  PLT  for  different  mean  processing  time, 
processing  time  causes  a larger  PLT  and  thus  more  holding  cost. 

2.  Increasing  penalty  cost.  With  increased  penalty  cost,  the  optimal  PLT  increases 
for  any  levels  of  holding  cost,  as  in  Figure  4.4  (a),  assuming  that  the  average  pro- 
cessing time  remains  a constant.  It  is  obvious  that  when  penalty  costs  are  high, 
it  is  inclined  to  have  a larger  PLT  to  avoid  a high  penalty  cost.  Furthermore,  we 
also  notice  that  the  optimal  expected  cost  is  proportional  to  the  penalty  cost  (see 
Figure  4.4  (b)).  In  other  words,  high  tardiness  costs  cause  high  optimal  expected 
costs.  This  is  true  because,  given  that  any  other  parameters  are  unchanged,  high 
penalty  or  holding  costs  result  in  high  overall  costs. 

3.  Increasing  holding  cost.  In  Figure  4.5  (a),  note  that  as  the  holding  cost  is 
increased,  the  optimal  PLT  decreases  for  any  levels  of  penalty  cost.  One  explanation 
is  the  fact  that,  since  a holding  cost  occurs  when  the  PLT  of  an  order  exceeds  the 
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(a)  (b) 


Figure  4.4:  Optimal  PLT  and  overall  cost  for  different  penalty  cost,  (a)  PLT  vs 
penalty  cost,  (b)  Overall  cost  vs  penalty  cost.  A = 27. 

actual  job-completion  time,  the  optimal  PLT  is  forced  to  be  decreased  to  avoid  a high 
holding  cost,  given  that  other  parameters  are  remained  unchanged.  On  the  other 
hands,  as  we  can  see  in  Figure  4.5  (b),  the  optimal  cost  increases  as  the  holding  cost 
increases,  given  that  other  parameters  are  remained  the  same. 


4.6  Conclusions 

Stochastic  production  systems  with  limited  resources  are  discussed.  We  have 
modeled  the  problem  of  setting  optimal  planned  order  leadtime  in  a single  stage, 
make-to-order  production  system.  As  a solution  approach,  we  derived  the  distri- 
bution of  actual  process  completion  time  of  individual  orders  and  then  determined 
the  optimal  planned  leadtime.  We  also  identify  that  our  model  is  equivalent  to  a 
simple  newsboy  problem  when  there  is  only  one  customer  order  (the  newly  arrived 
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Figure  4.5:  Optimal  PLT  and  overall  cost  for  different  holding  cost,  (a)  PLT  vs 
holding  cost,  (b)  Overall  cost  vs  holding  cost.  A = 27. 

one  itself)  in  the  system.  For  the  case  that  multiple  orders  reside  in  the  system, 
we  have  shown,  through  the  computational  studies,  that  the  optimal  PLT  is  propor- 
tional to  a unit  penalty  cost,  but  is  reciprocal  to  a unit  holding  cost,  whereas  the 
optimal  expected  cost  for  our  lead-time  model  is  proportional  to  both  penalty  and 
holding  costs.  These  results  are  rather  acceptable  because  when  a tardiness  cost  is 
much  higher  than  a holding  cost,  manufacturers  try  to  offer  longer  PLT’s  to  their 
customers  to  avoid  expensive  penalty  costs,  and  vice  versa. 


CHAPTER  5 

CONCLUSIONS  AND  FUTURE  RESEARCH 

This  dissertation  investigates  issues  on  stochastic  systems  with  limited  resources. 
We  study  stochastic  multi-item  inventory  models  with  limited  warehouse  capacity 
for  the  cases  of  both  equal  and  unequal  replenishment  intervals.  We  also  discuss 
stochastic  production  systems  in  which  a single  machine  serves  a series  of  customer 
orders.  The  primary  contributions  of  this  dissertation  are 

• The  development  of  heuristics  to  estimate  replenishment  quantities  when  a 
warehouse  confronts  a limited  capacity  and  individual  products  have  unequal 
replenishment  intervals. 

• The  derivation  of  the  optimal  order  leadtime  in  a single-stage,  make-to-order 
production  system. 

We  first  characterize  a stochastic  multi-item  inventory  model  for  the  case  of  equal 
replenishment  intervals  with  limited  warehouse  capacity  by  using  dynamic  program- 
ming. We  use  induction  to  derive  the  convexity  of  objective  cost  function  and  show 
that  a myopic  replenishment  policy  is  optimal  for  both  a finite  and  an  infinite  plan- 
ning horizon.  Furthermore,  the  optimal  policy  in  each  period  is  easy  to  calculate 
due  to  the  independence  of  the  one  period  expected  cost  over  periods. 
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For  the  unequal  replenishment  interval  case,  we  completely  characterized  the 
optimal  replenishment  policy  for  the  two-item,  two-period  case.  For  more  general 
cases,  because  of  the  inability  to  show  convexity  of  the  objective  function,  we  con- 
struct three  heuristics  to  determine  replenishment  quantities  for  the  case  of  unequal 
replenishment  intervals. 

Numerical  testing  of  these  heuristics  suggests  that  they  yield  near  optimal  so- 
lutions for  a small  set  of  problems.  The  performance  of  these  heuristics  has  been 
evaluated  for  a number  of  problems  with  various  cost  and  demand  parameters.  When 
there  is  enough  slack  in  the  resource  constraint,  all  heuristics  demonstrate  identical 
results. 

The  first  heuristic,  denoted  Heuristic  A,  which  uses  a separate  capacity  for  each 
product  for  replenishment,  is  simple  enough  for  practical  use,  but  is  not  able  to  detect 
a scarce  resource  well,  compared  to  the  other  two  heuristics,  in  most  problems.  This 
behavior  is  plausible  because  the  heuristic  is  more  likely  to  retain  unused  capacities 
for  some  products  while  some  others  suffer  from  the  lack  of  the  resource. 

Heuristic  C adopts  a replenishment  policy  which  replenishes  products  up  to 
the  total  available  capacity  of  the  system  at  the  time  instant  of  replenishment. 
This  heuristic  outperforms  Heuristic  A in  most  cases  because  it  does  not  carry  the 
capacity-inefficiency  confronted  by  the  Heuristic  A , but  leads  to  poorer  performance 
for  the  problems  with  high  penalty  costs.  This  is  because,  in  this  situation,  Heuristic 
C is  too  aggressive  to  replenish  products  in  one  period  so  that  the  warehouse  suffers 
from  the  lack  of  resource  in  the  following  replenishment  period,  causing  high  penalty 


costs  to  the  warehouse. 
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On  the  other  hand,  the  second  heuristic,  denoted  Heuristic  B,  adjusts  separate 
capacities  in  each  replenishment  period,  based  on  the  availability  of  any  unused 
resource.  This  characteristic  discriminates  Heuristic  B from  Heuristic  A and  dimin- 
ishes the  capacity-inefficiency  of  Heuristic  A.  Furthermore,  the  capacity-separation 
scheme  of  the  heuristic  leads  to  better  performance  than  Heuristic  C. 

Stochastic  production  systems  with  limited  resources  are  also  discussed.  We  have 
modeled  the  problem  of  setting  optimal  planned  order  leadtime  in  a single  stage, 
make-to-order  production  system.  As  a solution  approach,  we  derived  the  distri- 
bution of  actual  process  completion  time  of  individual  orders  and  then  determined 
the  optimal  planned  leadtime.  We  also  identify  that  our  model  is  equivalent  to  a 
simple  newsboy  problem  when  there  is  only  one  customer  order  (the  newly  arrived 
one  itself)  in  the  system.  For  the  case  that  multiple  orders  reside  in  the  system, 
we  have  shown,  through  the  computational  studies,  that  the  optimal  PLT  is  propor- 
tional to  a unit  penalty  cost,  but  is  reciprocal  to  a unit  holding  cost,  whereas  the 
optimal  expected  cost  for  our  lead-time  model  is  proportional  to  both  penalty  and 
holding  costs.  These  results  are  rather  acceptable  because  when  a tardiness  cost  is 
much  higher  than  a holding  cost,  manufacturers  try  to  offer  longer  PLT’s  to  their 
customers  to  avoid  expensive  penalty  costs,  and  vice  versa. 

We  conclude  this  dissertation  by  pointing  out  possible  extensions  of  the  stochastic 
models  considered  in  this  study.  The  current  optimization  model  is  established  for 
the  fixed  warehouse  capacity  over  time.  One  interesting  development  is  to  consider 
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a variable  warehouse  capacity.  This  extension  is  important,  since  in  many  contexts 
the  warehouse  capacity  would  be  variable.  For  example,  extra  warehouse (s)  can  be 
leased  to  accommodate  more  inventories  to  satisfy  customer  demands.  Or  Federal 
Express  uses  its  giant  warehouse  for  product  replenishment  as  well  as  for  its  daily 
operation,  such  as  mail  classification.  The  total  available  warehouse  capacity  then 
varies  according  to  the  size  and/or  the  duration  of  daily  operation.  In  this  case, 
the  warehouse  capacity  in  period  j,  say  Vj,  varies  over  periods  and  we  redefine  the 
problem  in  (3.17)  to  determine  the  separate  capacity  of  each  item  for  Heuristic  A as 
following 


Find 

s.t. 


p*  = min {/ij  > 0} 


£ < Vj 

i= 1 


V-  ■ = Frl 

vh3  ri 


( bjpj  - Hj  \ 
\bi(hi  +Pi)J  ’ 


i — 1,2,  •••,  l 


where  Vy  is  the  separate  capacity  of  item  i in  period  j.  We  also  can  redefine  the 
processes  for  Heuristics  B and  C using  the  same  manner  as  above  to  determine  the 
replenishment  policies. 

Another  research  effort  may  continue  along  by  relaxing  the  assumption  of  zero 
delivery  leadtimes  between  manufacturers  and  a warehouse.  It  is  also  of  interest  to 
analyze  the  warehousing  problem  with  the  additional  assumption  of  nonzero  fixed 


ordering  costs. 


APPENDIX  A 


There  are  two  cases  to  solve  equations  (b)-(d)  and  (4)  in  the  single-item  case. 
Assuming  that  F_1  exists, 


(i)  : I + Q -Vo  ^ 0. 

In  this  case,  we  get  p = 0 from  (c),  and  then  the  following  is  immediate  from 
(4): 


(A.1) 


Combining  (b)  and  the  condition  (i),  we  have: 


I + Q-V0<  0 


or  F-1 


P 


h + p 


Vo  < 0 from  (A.l) 


(ii)  : I + Q — \ o = 0.  In  this  case,  it  is  obvious  that  Q = V0  — I (*)  from  the 
condition  (ii).  Then,  from  (2.4)  and  (c), 


p = p—  (h  + p)F(I  + Q)  > 0 


or  Vq  < F 


-l 


P 

h + p 


from  (*) 
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=»  Q = Vo -I 


if  Vo  <F' 


i / P 


h + p 


Combining  those  two  results,  we  have: 


Q*  = 


F~l  (*)  - 1 
V0 -I 


if  F-'(^)>V0 


Or  simply 


Q*  = min 


F- 


i 1 P 
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-I 


APPENDIX  B 


The  first  and  second  partial  derivatives  of  gj  with  respect  to  Qij  can  be  deter- 
mined as  follows: 
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(B.2) 


where  Ajj  = Iij  + Qij  and  Fij  is  the  cumulative  distribution  function  of  demand 
Dij.  Since  the  right-hand  side  of  (B.l)is  always  nonnegative  and  (B.2)  is  equal  to 
zero,  the  Hessian  matrix  Bj  is  positive  semi-definite. 
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APPENDIX  C 


We  use  the  induction  method.  First,  consider  the  period  m.  Define 
Lm  = Gm{I\,mi  ...,  Il,mi  Ql,m>  •••)  Ql,m)  "P  /fin  "P  Sj— 1 Qi,m  ~ ^0^  • Then, 

from  the  first-order  condition,  we  get 
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Solving  (C.l)  and  the  KKT  condition:  /im(X)-=1  + Z)!=i  Qi,m  - Vo)  = 0,  we  get 
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